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Abstract 

Short overview of experiments with SND detector at VEPP-2M e + e~ collider in 
the energy range 2Eq = 400 1400 MeV and preliminary results of data analysis 
are presented. 
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1 Introduction 

The SND detector was proposed in 1987 [@ to continue successful series of experiments 
with the ND detector at VEPP-2M e+e~ collider |§ in Novosibirsk. The goal of 
experiments was a detailed investigation of e + e~ annihilation processes in the same energy 
range 2E = 400 4- 1400 MeV, in particular the final states with photons and other neutral 
particles (ir°, T], u, K s ), decaying into photons. 

The main part of the SND is a spherical electromagnetic calorimeter, consisted of 
1630 NaI(Tl) crystals. The total weight of NaI(Tl) is 3.6 t, the solid angle coverage is 
~ 90% of Air steradian. Energy resolution for photons is &e/E = 4.2%/ t/E(GeV) 0, 
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angular resolution is about 1.5 degrees. The angles of charged particles are measured 
by two cylindrical drift chambers covering 95% of An steradian solid angle. Angular 
accuracy of charged tracks measurements (er) is about 0.4 and 2.0 degrees in azimuth in 
polar directions respectively. From the outside SND detector is covered by muon system, 
consisting of streamer tubes |]] and plastic scintillation counters. 

2 PHI-96 Experiment 

Detailed description of PHI-96 experiment and first physical results, based on about half of 
the statistics, were published in flU . The PHI-96 experiment was carried out in the period 
from February 1996 until January 1997. Seven successive scans PHL9601 -r- PHL9606 and 
PHL9608 were performed in the center of mass energy range 2Eq from 980 to 1044 MeV. 
Data were recorded at 14 different values of the beam energy. The total of 4.4 pb~ l 
integrated luminosity was collected, corresponding to 8.3 million mesons produced. 
Most of the results presented in this paper are based on full recorded statistics. 

3 MHAD-97 Experiment 

The MHAD-97 experiment was described in 0. In the period from January until June 
1997 two scans were performed in the energy range 2Eq from 960 to 1380 MeV with a 
step of 10 MeV and total integrated luminosity of 6.3 pb^ 1 . At present the MHAD9701 
scan is accessible for processing completely, while only part of MHAD9702 scan has been 
preprocessed. The results of this paper are based on the integrated luminosity 2.8 pb~ l 
from both scans, corresponding to production of about 1.6 ■ 10 5 p + p~ pairs. 

4 PHI-97 and OMEGA-98 Experiments 

The 1997 — 1998 data taking runs started in October 1997. The goal was a significant 
increase of the meson statistics in order to improve the accuracy of measurements of 
electrical dipole decays — > 7r°7r°7 and — > i]7T 0, y, observed with SND for the first time. 
Three big scans PHL9703, PHL9801, PHL9802 were carried out, covering 16 energy 
points in the center of mass energy range 2E = 984 -h 1060 MeV. The total recorded 
integrated luminosity was 10.3 pb' 1 , corresponding to 15 million mesons produced, that 
is two times more than in PHI-96 experiment. 

From March to June 1998 both SND and CMD-2 detectors took data in the energy 
range 2E from 970 MeV down to 360 MeV (OME_9803 scan). The goal of this run was 
measurement of hadronic production cross sections and search for rare decays of p and 
uj mesons. Integrated luminosity of 3.4 pb~ l was collected, corresponding to ~ 4.4 ■ 10 5 
p + p~ pairs and ~ 1.9 • 10 6 p mesons produced. Significant part of data — 1.2 were 
accumulated in the energy range 2E = 760 4- 800 MeV, where ~ 1.2 • 10 6 to mesons were 
produced. In Fig. [I] a week-by-week graph of statistics, recorded by SND is shown. Fig. ^| 
shows the VEPP-2M luminosity averaged over experimental time as a function of energy. 
Solid curve in Fig. [| corresponds to Eq dependence, normalized to the experimental value 
L = 2 ■ 10 30 cm" 2 s _1 at E = 500 MeV. All listed above experiments were carried out 
with a superconducting wiggler [|J. 
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A distribution of full integrated luminosity accumulated with SND up to now is shown 
in Fig. p. The data processing of last runs PHI-97 and OMEGA-98 will be started in the 
end of 1998. We believe, that preprocessing of this data will be finished in half a year. 
After that, new data will be available for analysis. 

5 General remarks on data processing 
5.1 The analysis of "neutral" modes 

SND has an advantage in detection of pure "neutral" final states. Up to now the PHI- 
96 and MHAD-97 experimental data were analyzed in order to study or search for the 
following 0-meson decays: 

e + e- ^(f)^ 7r°7r°7 -> 57, (1) 

e+e~ — > 4> — > ?77r 7 — > 57, (2) 

e+e~ — ► (f) — > ?77r 7r 7 — > 77, (3) 

e+e" -> -> r]/(958)7 -> ^7r°7r°7 (4) 

e+e~ — > — * 777 — * 37r°7 — > 77, (5) 

e+e~ — > — > KsKl — > neutrals, (6) 

e+e~ — > — > 777 — > 37, (7) 

e+e" -> -> 7T°7 -> 37, (8) 
as well as non-resonant or non 0-meson processes: 

e+e" -> CU7T° -> 7T 7T 7, (9) 

e+e" — > 77, (10) 

e + e-^37. (11) 

e+e" ^4 7 , 57. (12) 

e+ e - -> K S K L , (13) 

e+ e - -> K 5 iT L7 , (14) 

e+e" - KT^tt , (15) 

e+e - — > 777 — > 37T°7, (16) 

e+e" -> if 5 K L , K 5 -> ttV, (17) 

e+e" -> cjttV -> 3tt°7, (18) 

e + e- -> / 2 (1270) -> ttV. (19) 

In analysis of any processes listed above, all remaining processes should be considered 
as a background. Processes with any number of photons in the final state could mimic any 



certain process, particularly in case, when the cross-section of the background process is 
much larger. The number of found photons in an event could be less than the number of 
produced photons, due to the merging of close showers in the calorimeter, or when some 
of particles escape detection. On the other hand, additional photons could be found in 
events due to shower splitting or beam background photons hitting the detector. 

Parameters, which are widely used in different analyses are listed below: 
iV 7 — number of found photons. 
Etot — total energy deposition. 
Ptot — absolute value of total momentum. 

E tot /2E Q — total energy deposition in the calorimeter, normalized by the center of mass 
energy. 



Ptot I 'E to t — absolute value of total momentum of all detected particles in assumption, 
that all particles are electrons and photons, normalized by the total energy deposition. 
E np /2E — total energy deposited by neutral particles, normalized by the center of mass 
energy 2E . 

Ei — energy deposition of z-th particle. 

E~f max — maximum energy of neutral particle in an event. 

9i — polar angle of z-th particle (particles are enumerated in the following way: charged 
particles first, then neutral particles in descending order in energy within each group). 
Omin — minimum polar angle between particle and beam direction in an event. 
4>i — azimuth angle of z-th particle. 

( — the likelihood of a hypothesis, that given transverse energy profile of a cluster of hit 
crystals in the calorimeter can be attributed to a single photon [fjH , ||10|| . This parameter 



allows one to separate events with isolated photon showers, from events, which have 
overlapping showers or group of hit crystals from i^-meson decay or nuclear interaction. 
X% — parameter, describing the degree of energy-momentum balance in an event under 
assumption that all particles are photons and electrons. 

xj — parameter, similar to x\i describing likelihood of assumption, that there is an 
intermediate state / in event (for instance / could be 7r°7r°7, riii ^, etc.). 
N w o — number of found 7r°-s. 

For example for primary selection of 5-7 events of the processes (^]), (0), and (0), the 
following cuts were imposed: 



AT 7 = 5; N cp = 0; 
0.8 < E tot /2E < 1.1; P tot /E tot < 0.15; (20) 

@min ^ 27 . 

5.2 Analysis of the processes with charged particles 

Up to now the following processes with charged particles in the final state were analyzed: 

e + e~ — > (j) — > i]e + e~ — > e + e^ r yy, (21) 

e + e~ — > 4> — > r?7, r\ — > e + e~7, (22) 

e + e~ — > (f> — > rj'j, r{ — > n + n^r], r\ — > 77, (23) 

e + e~ — > — > 7r + 7r~7r + 7r~, (24) 
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e+e~ -> wtt° -> tt+tTttV, (25) 

e+e - — > e + e~7; (26) 

e+e - — > e + e~77; (27) 

e+e" -> -> (28) 

e+e" - -> i^ L , Ks^tt+tt-; (29) 

e+e" -> tt+tt-tt ; (30) 

e+e" -> -> tt+tt-tt ; (31) 

e+e^ — > — > 777, 77 — > 7r+7r-7r , (32) 

e+e" -> 7r+7r-7r°vr . (33) 

e+e - — > 7r+7r-7r+7r- (34) 



Let us list typical parameters, used in analyses of events with charged particles: 
N cp — number of detected charged particles; 

Ri — distance between z-th particle track and beam axis in R — plane; 

Zi — Z coordinate of the point on the track of z-th particle, closest to the beam axis; 

N W ire — number of fired wires in the drift chamber. 

a,ij — spatial angle between z-th and j-th particles. 

dE/dxi — energy deposition in drift chamber for z-th particle. 

Parameters of kinematic fit like Xe are a ^ so used with additional assumptions about 
masses of charged particles in a similar way as x% 



6 Physical results from the PHI-96 experiment 

6.1 The (f) — > 777 decay 

The — > rj'-f decay is a radiative magnetic dipole transition of into rj meson, studied 



previously in many experiments [IT]. In this work we measured the — >• 777 decay rate 
into multi-photon final state: e+e~ — > — > 777 — > 37r°7 — > 77 (|[). To suppress the 
background the events were selected satisfying the following criteria: 

• iV 7 = 6, 7, 8; N cp = 0; 6 mm > 27°; 

• P/2E Q < 0.15; 0.8 < E tot /2E < 1.1; X \ < 30. 

The presence of tt° mesons in event was not required in the kinematic fit. As a result of 
such selection cosmic background was completely rejected and main background process 
(0) was suppressed significantly, as one can see in the distribution in m rec . 7 , recoil mass 
of the most energetic photon in the event (Fig.|27|), where peak at 77 mass dominates. For 
final selection of 777 events the soft cut on m rec . 7 was used: 400 MeV < m rec . 7 < 620 MeV. 
The estimated detection efficiency is close to 10% and background contribution, obtained 
using the number of events from the m rec . 7 interval 620 4- 840 MeV, is about 2%. The 
ratio of the background events numbers in these two regions was taken from simulation 
of the process @. The fit of visible cross section was done for each of the 6 scans 
separately. Free parameters of the fit were the shift of the energy scale of VEPP-2M, 
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meson width, and the branching ratio B(<f) — > 777). Also taken into account were collider 
beam energy spread, instability of average beam energy, the radiative corrections, the 
VDM contribution of p and u resonances, and presence of dead calorimeter channels. 
The fit results are listed in Table |l|. Averaging the data from the Table |l] one can obtain 

Table 1: The fit results for the <fi —>■ 777 decay in 6 independent scans. 



Experiment 


^events 


IV MeV 


5(0^ 777),% 


PHL9601 


1045 


4.34 ±0.43 


1.141 ±0.066 


PHL9602 


1436 


4.00 ±0.34 


1.188 ±0.044 


PHL9603 


2163 


4.29 ±0.33 


1.192 ±0.036 


PHL9604 


1241 


4.53 ±0.39 


1.127 ±0.059 


PHL9605 


2222 


4.38 ±0.24 


1.271 ±0.049 


PHL9606 


1709 


3.97 ±0.27 


1.339 ±0.059 



the branching ratio 

B((f) -> 777) = (1.209 ± 0.028 ± 0.050)%. 

Here the first error is statistical and the second is a systematic one, estimated to be 4.2 %. 
The systematic error is mainly determined by the following contributions: 

• the systematic uncertainty of normalisation (3%); 

• the background subtraction error (1%); 

• the error in the detection efficiency (1%); 

• the error in B(<p -> e + e~) (2% fly]); 

• the error in B(r] -> 3vr°) (1.2% |[TJ); 

• the error in value of p — (f> interference term (1%). 

The results in Table [1] show some difference between separate PHI-96 scans. The scale 
factor for the branching ratio, calculated according to PDG recommendations, is equal to 
1.4. It was taken into account in the presented statistical error. 

The result obtained in the present work is in agreement with PDG value (1.26±0.06)% 
IT]. At the moment it is the most accurate single measurement of B{<p — > 777). 



6.2 Analysis of the — > r]e + e — ► 77e + e decay 

The decay 4> — > r/e + e~ is closely related to the radiative decay (f> — > 777, where, instead of 
a real photon, virtual one is produced, decaying via the channel 7* — > e + e~ . This type of 
decays is called Dalitz or conversion one. Experimental study of these processes is a test 
of quantum electrodynamics. Besides this, existence of transition form factor [T2|] affects 



invariant mass spectrum m e + e -, which was evaluated in different models. Theoretical 
value of B((p —>■ r]e + e~) for a unit form factor is 1.1 ■ 10~ 4 



This decay was observed for the first time in our experiment with the detector ND [[13 
in 1985, with a branching ratio of (1.3i ;g) • 10~ 4 . Later similar result was obtained with 
the CMD-2 |1 detector: B((f) -> rje+e~) = (1.10 ± 0.49 ± 0.19) • 10" 4 . 



In present analysis the — > 7/e + e — > 77e + e decay was studied in the reaction 
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e^e 



r]e + e — > e + e 77 The following event selection criteria were used: 



iV 7 = 2, N cp = 2, 

E tot /2E > 0.8, R h R 2 < 0.2 cm, 
10° < ai )2 < 110°, L liee < 6, 

where L liee is a logarithmic likelihood function, obtained as a result of kinematic fit, 
using the energies and angles of all particles (invariant mass of two photons m 77 being 
evaluated as well). All 7 0-meson scans were processed with these selection criteria. 

The distributions of experimental events over invariant mass of two photons together 
with corresponding distributions of simulated — > r]e + e~ and QED e + e~ — > e + e~77 
events are shown in Fig.[|. The simulated distributions are normalized to the branching 
ratio B((p — > i]e + e~) = 1.3 • 10~ 4 from [111] for — > r]e + e~ , and to the total integrated 



luminosity for QED e + e~ — > e + e~77 reaction. One can see that the experimental back- 
ground is almost entirely provided by QED process, which has identical final state as the 
one, searched for. Now let us evaluate the result in two ways. 

First, the experimental distribution over m 77 was fitted with a background, approx- 
imated by a third order polynomial and Gaussian, representing the — > r/e + e~ process. 
As a result of the fit (Fig.|5|) we obtain the number of events in the peak Nq = 54 ± 15. 
Branching ratio reads: 

B{6 -> 7?e + e-) = -i = (1.42 ± 0.39) • 10" 4 , (35) 

where = 8.3 • 10 6 is the number of produced 0-mesons, e = 4.45 ± 0.18% is a detection 
efficiency, obtained by simulation. Systematic error in e here was neglected. 

Second, let us derive branching ratio — > r/e + e~, using the energy dependence of 
visible cross section of the events, satisfying the additional condition 500 MeV < m 77 < 
600 MeV. The fitted resonance curve is shown in Fig|| 0-meson mass and width were 
fixed at the values = 1019.41 MeV and 1^ = 4.43 MeV [[TTJ] . As a result of the fit 
the following parameters values were obtained: B(<f) — > 7/e + e~) = (1.71 ± 0.42) ■ 10~ 4 , 
background level cr& fl = 8. 7*2.7 pb- Apparently, the value of B{<p — ► r]e + e~) is excessive 
due to resonant background admixture (supposedly a greater part is due to — » 7r + 7r _ 7r ). 
Correction for this resonant background we estimate in the following way. Let £ be the 
ratio of the resonant and non-resonant background level. The values of £, derived by 
fitting a resonance curve to energy distributions of background events of the intervals 
m 77 = 350 ± 50, 450 ± 50, 650 ± 50 and 750 ± 50 MeV, are well approximated with a 
straight line. Thus at the value of m 77 ~ 550 MeV we obtain 

£L 77=55 o = (5-8±2.7)-10- 2 . 

After subtraction of the resonant background one gets: 

B((p^r]e + e-) = (1.21 ±0.51)- 10~ 4 (36) 

Obtained in essentially different ways the two values of B(<f) — > r]e + e~)\ (35) and (36) 
are in a good agreement. As a final result one can choose (35) as the most precise of them 
and consider their difference as a measure of systematic error of ±0.21 or ±15%. After 
taking into account other known sources of systematics the error estimate increases up to 
16% and the final result reads as 

B(<f> -> rie + e~) = (1.42 ± 0.39 ± 0.23) • 10" 4 . (37) 

It is close to the theoretical value and does not contradict previous measurements. 
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6.3 Study of the process 



777 , 77 — ► e T e 7 



The Dalitz decay 77 — ► e + e 7 was studied in the reaction e + e (f> — ► 777, 77 — > e + e 7 (|22|). 
The ratio of probabilities of the Dalitz decay and the two photon decay of 77 meson is 



given by the following expression] 12 



d B(r) — ► e + e 7) 
(im e 2 e 5(77 — > 77) 



2a |F,K e ) 



3tt 



m; 



2m? 



771; 



1 - 



4m; 



1/2 



m; 



m; 



mf 



- 77-meson transition form factor. The main 
? J ~ 1 . The expected branching ratio is: 



where m ee is a e + e invariant mass 
contribution is given by small m ee , where 
B(r] -> e+e"7) = 6.3 • 10" 3 . 

This decay was measured in only one experiment [37], where in two successive papers 
the results differ by a factor of 3. We believe, that our measurement could clarify the 
situation with this decay. 
Data analysis 

The following event selection criteria were used : 
• positive charged event trigger; 



N, 



cp 



2; 2 < AL < 3; R h R 2 < 0.5 cm; \Z X \, \Z 2 \ < 10 cm; 



• 6 min > 36°; A0 ee > 5°; E tot /2E > 0.8, P/E tot < 0.15; X % < 15; 

• the events with invariant mass of photon pairs close to masses of 7r° 77 mesons 110 MeV < 
m 77 < 170 MeV; 500 MeV < m 77 < 600 MeV; are excluded. 

• E 1 m i n > 50 MeV; recoil mass of one photon is close to 77 meson mass in the decay 
^ 777 : 0.6 < E 7 /E < 0.7 ; 0.1 < m e+e -/m ri < 0.7. 

After application of these selection criteria 60 events survived. We estimated their 
origin in the following way: 



Experiment (total) 

Contribution from the process e + e~ — > e + e~77 
Contribution from the process — > r]e + e^ 
Contribution of the process under study 



60 
35 
0.5 
24.5 



Taking the efficiency 4.15% from simulation we have 

24.5 



e e 



'l) 



0.0415 • N$ ■ B(4> -> 7/7) 



6.1 • 10" 



The energy dependence was fitted by the sum of QED non-resonant background and 
resonant contribution from the decay under study (fig.0). The following fit result was 
obtained B(rj —>■ e + e~j) = (6.8 ± 1.1) • 10~ 3 . The systematic error was estimated to be 
about 10%. The final result is: 

5(7/ -> e+e^) = (6.8 ± 1.1 ± 0.7) ■ 10~ 3 . 

The result obtained in this work is in good agreement with the PDG value (4.9 ±1.1) • 10~ 3 
| n|] and with the prediction 6.3 • 10~ 3 . The accuracy of our measurement is close to the 
world average value |[T] . 
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6.4 Search for the — ► r]7T 7T j decay 

Radiative decays of vector mesons V —>■ Mj, where M is a scalar or pseudo-scalar state, is 
an important source of information on the structure of these states. For V = 0(1020) only 
main decays — > 777, 7r°7 [jn| were studied by now at a relatively high level of accuracy 
of < 10%. Regarding rarer decays like — > t/(958)7, /o(980)7, a (980)7, even the latest 
data |16|, |10], |17| are still preliminary and have rather low accuracy. 



The — > 7/7 decay was observed in the only experiment with CMD-2 detector at 
VEPP-2M collider [T^J with a branching ratio of (1.2^5) • 10 -4 . The quark model pre- 



diction for this value |l8j is (0.7 -4- 1.0) • 10 -4 under assumption that there is no gluonium 
component in rj' . In case of rj' consisting of pure gluonium, the expected branching ratio 
of the decay — > 7/7 would be 4-10~ 6 |i~9|. 

In this work the decay — > 7/7r 7r°7 — > 77 was studied. The pseudo-scalar M = 777r 7r° 
final state could appear from the decay of 77' (958) or from states with higher masses 
and large widths like 7/(1295) or 7/(1440). The expectation for the — ► 7/7 — > 777r 7r°7 
branching ratio, based on quark model predictions is B{<p —>■ 7/7) ■ B(rf — > t]tt tt°) = 
(1.5 -=-2.0) ■ 10" 5 . 

The selection criteria In this work the reactions (0) and (ffi) were studied, both 
producing a 7-photon final states. The main background comes from the 0-meson decays 
© and (§. 

The events of the processes under study @ and (§) were selected using the following 
criteria: 

Sel.l - event configuration does not contradict energy and momentum conservation: 
Xe < 15. This requirement suppresses the background from the process (f|). 

Sel.2 - the energy of the most energetic photon in the event E max /E < 0.65. Fig.[| 
shows, that such a threshold almost completely rejects the background from process (|5]), 
although the detection efficiency for the process (|j) becomes twice smaller. 

Sel.3 - the transverse spread of each shower in the calorimeter does not exceed certain 
standard value, characteristic for individual showers ( < 0. This cut suppresses the 
background events from the process (^), where some close photon showers merge together. 

Sel.4 - among all photon pairs in the event, three meson candidates must be found: 
two 7r°-mesons with the effective mass 110 < m 77 < 160 MeV and one 77-meson with 
520 < m 77 < 580 MeV. Analysis of simulated events of process (U) showed, that the 
pair of most energetic photons in the event is produced by 7/(550) — > 27 decay. This fact 
significantly reduces combinatorial background in the search of tt° and 77 candidates in 
the event. 

Sel.5 - the effective mass of T/7r 7r°-system is close to rj' mass. The obtained mass 
spectra of rj' for simulated events of process (^) is shown in Fig|]. 

The data analysis. The — > 7/7 decay detection efficiency with respect to all events 
of this decay e\ = (0.73 ± 0.09)% was obtained using simulation of the process (Q) with 
the selection criteria Sel.l — Sel.5. Corresponding efficiency with respect to events of the 
process (H) with exactly 7-photons in the final state is much higher ~ 9%. Among all 
experimental events N x =3 were found, satisfying selection criteria. They are shown as 



black triangles in the scatter plot Fig.[LC| The vertical axis of the plot depicts the effective 
mass of the found 7/(550) meson, while the horizontal axis shows its momentum. After 
the additional cut on 7/(550) momentum P v > 50 MeV ' jc two experimental events survive. 
Simulated events of the process (|j) are shown in the Fig.|TD| as circles. One event, shown 
as a square, was found in the analysis of the simulated — > 7/7 background sample. The 
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estimated detection efficiency for this background is close to 10~ 5 . 

In the analysis of 4.7- 10 5 simulated events of the reaction (§) no events of <p — > K$Kl 
decays, satisfying selection criteria Sel.l — Sel.5 were found, but one should take into 
account, that <p ~ > K$Kl simulation statistics is 5 times smaller than the experimental 
one. Because the simulation is not precise, certain contribution of the process ([]) into the 
sample of two selected experimental events cannot be completely excluded. 

Both found experimental events correspond to a collider beam energy near 0(1020) 
maximum, which confirms their resonant origin. If both events originate from the process 
(H), the branching ratio would be equal to (0.4 ± 0.3) • 10 -4 . But as it was pointed out 
above, the background origin of these events from processes (||) and (Q) could not be 
excluded. Thus only upper limit on the branching ratio of the decay (f|) could be reliably 
established. At a confidence level of 90% it is equal to 

Br (<p - rj'j) < = 1.1 . lO" 4 , (38) 
where coefficient k = 2.7 comes from Poisson distribution of observed decays [|TT |. 



To check the whole procedure on possible systematic errors we studied another decay 
chain ([|) with a similar 7-7 final state. Its probability is known and could be compared 
with an obtained value. The selection criteria Sel.l<S>Sel.3 were imposed. In addition it 
was required for each selected event to contain three 7r° candidates. We found B((f> — > 
777) = 1.05 • 10~ 2 , that is 18% lower than the table value. Such large discrepancy could 
be attributed to excessively stringent selection criteria Sel.l®Sel.3, used in the search 
of process (§). Now we can correct the detection efficiency e\ obtained using simulated 
events by 18% down to a value of e\ = 0.6%, and use it within the expression (^) as an 
experimentally obtained detection efficiency. 

In search for the process (|3|) without any assumptions about its dynamics the cuts 
Sel.l — Sel.4 were applied. The detection efficiency £2 was determined from a sample of 
simulated events of the reaction (|3])with the matrix element similar to that of the process 
(d) and effective masses of 777r°7r° system ranging from 850 to 1000 MeV. The obtained 
value of €2 is nearly constant in this effective mass range and equals to ~3% with the 
correction described above imposed. In the experimental sample N x =3 events were found, 
one event was found in simulated sample of background process (|5|) and none were found 
in the simulated sample of (||). After formal background subtraction the estimated the 
number of events from the process @ is N x = (3 — 1) ± a/3 + 1 = 2 ± 2. Given such a 
small number of events and not reliable estimation of background, only upper limit of the 
branching ratio can be placed using the expression (^8|): 

B(<p -> 777r Q 7r°7) < 2 ■ 10" 5 . 

Conclusion. In the analysis of 7-7 events the following upper limits at 90% confidence 
level were placed: 

5(0 -+7/7) < 1.1- 10~* 

B{(j> -> 777r°7r°7) < 2 • 10~ 5 (incl.7/7). 

6.5 Observation of the — ► 7/7 decay 

The result of the first measurement of <fi — > r]' / ~f branching ratio in 7/ — > ?77r + 7r~, 77 — ^ 77 
decay mode with the CMD-2 detector at VEPP-2M [O is B(<p -> 7/7) = (1.2±^) • 10~ 4 . 



At the same time in SND experiment only upper limit was obtained [|7| in another decay 
mode: 7/ — > 7/7r + 7r _ , 77 — > 7r + 7r~7r°: B(cf) — > 7/7) < 1.7 • 10 -4 at 90% confidence level. 
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This work is devoted to the search for the — > 7/7 decay in the decay mode e + e — > 
rj /r y, rj' — > i]7r + iT~ , rj — > 77 (jffi). The background processes are the following: e + e~ — ► 777, 
77 -»■ tt+tt-tt fl32|); e+e" -> -> tt+tt-tt e+e" -> -»• cjtt -»• 7r + 7r-7r 7r fl2|). The 

main problem in search for (|23|) is, that the process (0) has the same final state, while its 
cross section is two orders of magnitude higher. Processes (|31"| , p5| ) give similar signatures 
and could contribute to the process under study. 

Event selection. Full statistics of 7 scans was analyzed, corresponding to the production 
of 8.3 • 10 6 (p mesons. The expected numbers of experimental events as well as numbers 
of simulated events are listed in the Table |2|. The following cuts were applied for the 



Table 2: The expected numbers of experimental events and numbers of simulated events 



for investigated and background processes. 



Process 


Expected number 
of experimental 
events — iV ex .p. 


Number of 

simulated 

events 

Nmnd. 




4> -> 7/7 


830, at 5(0 -> 7/7) = 10" 4 


9996 


0.083 


4> -> 777 


104580 


368176 


0.284 


(f) — ► 7r + 7r~7r° 


1.29 • 10 6 


75000 


17.3 


e + e~ — > ujit 


38 • 10 3 


33784 


1.12 


(j> -> K + K~ 


4.1 • 10 6 


24991 


163 


- K$Kl 


2.8 • 10 6 


194551 


14.5 



preliminary event selection: 

• charged particles trigger; 

• N cp = 2; iV 7 = 3; R\, R2 < 0.3 cm; \Zi\, \Z 2 \ < 6 cm; 8 < N wires < 12; 

• Xr < 50; xl < 50; 

• 6 min > 27°; < 145°; E tot /2E < 0.8; E NP tot /2E > 0.4. 

The cut on a n7T allows to exclude the background process e + e~ — > (f> — > K S K L , K s — > 
7T + 7r _ , — > neutral particles and possible background from Bhabha scattering and 
higher order QED processes. 21103 experimental events passed these cuts: 

The analysis of the kinematics of the reaction (|23|) shows, that the photon from the rj 
radiative decay always has the smallest energy E^ ~ 60 MeV, while the energy of photons 
from the decay rj — > 77 lies in the interval 170 MeV < E 7 i, E l2 < 440 MeV. In the pro- 
cess (|32"D the radiative photon has a maximum energy E yl ~ 362 MeV, while the energy 
of photons from the decay 7r° — > 77 lies in the interval 12 MeV < E l2 , E l3 < 364 MeV. 
It should be mentioned, that a probability for photons to reach the boarders of mentioned 
intervals is small due to the phase space suppression. Hence, there is almost no combi- 
natorial background for these processes, even taking into account the finite resolution of 
the detector. 

In Fig. |ll] the distribution of 11239 experimental events over recoil mass of the most 
energetic photon is shown. It was obtained with the additional cut \rn23 — Tn n o\ < 35 MeV 
(ir° meson mass resolution is equal to o m ~ 12 MeV). The clear peak on n mass from 
the process (^) is seen together with a broader distribution from the processes (|31~1 , |25|) . 
Fitting the data by the sum of Gaussian distribution and 7-th order polynomial one could 
estimate the number of events in the peak N V1 = 3623 ± 60. The order of polynomial 
is not important, if the mass range is limited to m rec7 < 700 MeV. The detection 
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efficiency is equal to e^Bij] — > n + n n°) = 3.63 ± 0.04%. It leads to the branching ratio 
B((p — > 777) = 1.20 ± 0.03%, where the error is a pure statistical one. The obtained 



result is only 5% smaller than the table value |Tl|, confirming validity of our analysis 
and absence of significant systematic errors. With additional cut, excluding 7r° meson 
|m 23 — m w o\ > 35 MeV, the peak at rj mass vanishes. So to exclude the events of the 
process (|32]) the cut 



^■ m ree .i-rra^ ^2 _|_ / TOag^m^g ^2 ^ ^ 

was imposed, where each denominator is equal to three times the resolution (3<r) for 
corresponding meson mass. 

In order to reach more effective suppression of background events, the characteristic 
kinematic features of the investigated process were used: the sum of three photons energy 
exceeds 607 MeV; the photon energy from the decay 77 — > 77 exceeds 170 MeV. For 
16512 experimental events the distribution over the energy deposition of neutral particles 
and over the reconstructed energy of the 2-nd photon is shown in Fig. |12| and Fig. [13] 
as well as estimated contribution of background processes. These figures show that the 
processes ( |3"T| ) and (E5[) are dominant. To reduce their contribution the following cuts 

• E NP tot /2E > 0.5, E l2 /E > 0.45, 

were applied, shown by arrows in Fig. O and Fig. [TBI. A cut ensuring the quality of 
photon showers in calorimeter was also used 

• C<o. 

367 events passed the cuts listed above. 

A scatter plot of invariant masses of photon pairs m 13 , m 2 3 and m 12 (for example see 
Fig. [TJ|) shows a concentration of events along the lines corresponding to tt° masses, which 
also confirms their background nature. To exclude events of this kind the cuts 

• |^i3 — m n o\ > 35 MeV, |m 2 3 — m^o | > 35 MeV 

were added. Second of these cuts simply tightens the condition, used for event selection 
of the decay (p — > 777. 

Scatter plot of 28 remaining experimental events in recoil mass of the photon with 
minimal energy — m rec ^ versus invariant mass of pair of most energetic photons — m\2 
is shown in Fig. |T5|d. The same distribution of simulated events for searched and main 
background processes is shown in Fig. |T5|a,b,c. Numbers of simulated events are not 
normalized. In order to display contributions from the processes (^T|) and fl2"o]), events in 
Fig. |15| c were selected without the cut on E l2 /E parameter. 

Projection of two dimensional plot on the axis 777,34 is plotted in Fig. |16|. Calculated 
contribution from the process (|23| ) at B(<p — > 7/7) = 10~ 4 , as well as the main background 
process Q3"2| ) are plotted in Fig. [TJ]. Besides the peak at 77 mass, background events of the 
processes fl3T|) ( p5|) form the peak at tt° mass. 

In Fig. [17| the second projection of two-dimensional scatter plot on the axis m rec ^ is 
shown with additional cut \rn12 — < 35 MeV. The peak at 7/ mass is observed, con- 
firming the existence of the <fi 7/7 decay. Histogram and smooth curve show distribution 
of 14 experimental events and its optimal fit by the sum of linear function and Gaussian 
with fixed parameters: = 957.5 MeV and a mrj , = 6.2 MeV, obtained from simulation 
of the process (|2"3"p. Hatched histogram and dashed line show the estimated contribution 
from the background process (|32|) and its approximation by the linear dependence. One 
can see, that the difference in the background estimation from experimental data and 
simulation does not exceed 0.6 events in the interval m v ± 20 MeV (2.6 and 3.2 events 
respectively). The number of useful events in the peak equals to iVy 7 = h.2~^%. The 
detection efficiency of the final state is equal to e^^ = 5.5 ± 0.6%. Taking into account 
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table values of B(rj' — > ?77r + 7r ) and B(r) —> 77) p| , we obtain 



5(0 - 7/7) = (6.71^) • 10- 5 , (39) 

which is nearly twice lower than CMD-2 result [117], but does not contradict it because of 
large error in both measurements. 



6.6 Study of the — > 7r tt 7 decay 

Search for — > 7r°7r°7 decay was first carried out in the ND experiment at the VEPP-2M 
e + e~ collider in 1987 and the upper limit B(0 — > 7r°7r°7) < 10~ 3 |2U], [T5|] was placed. As 



it was shown later by Achasov |^TJ , study of this decay can provide a unique information 



on the structure of the light scalar / (980) meson. Subsequent studies J22], |24|, [25], 



2T| , |28|| proved this idea. In these works different models of the / (980)-meson structure 



were considered. The most popular were two-quark model p9| , four-quark model 3(| 
and a molecular model [|25] . 

In 1997 the first indications of the process e + e~ — > — > 7r°7r°7 ([!]) were reported by 
SND []10| . In this work the results of analysis of full statistic of this experiment are 



presented. Main resonant background to the decay (fj) comes from the neutral decay 
e + e~ — > — > 777 — > 37r°7 (||) due to the merging of photons and/or loss of photons through 
the openings in the calorimeter. The main source of non-resonant background is a process 
e + e~ — > wit — > 7r°7r°7 @. The background from the — > pir — > 7r°7r°7 decay turned out 
to be small p4| , |31| , nevertheless, its amplitude was taken into account in the simulation 
of the process @. The background from the QED 5-7 annihilation process was found to 
be negligible. 

In order to suppress the contribution of background the events were selected satisfying 
the following criteria: 

• N y = 5; N cp = 0; 6 min > 27°; 

• 0.8 < E tot /2E < 1.1; P/2E < 0.15; C < 0; xl^ < 20. 

Parameter ( facilitates efficient separation of the events with well isolated photons from 
the events with merged photons or those produced by Kl mesons. The — > KgK^ 
tt 7t°Kl decay can contribute due to nuclear interactions of Kl mesons in the material of 
the calorimeter but after described cuts it is not seen at present level of statistics. The 
X^.o 7r o 7 , kinematic fit parameter, describes the degree of energy-momentum conservation 
in the event with additional requirement of presence of two tc° mesons. During this fit 
all possible combinations of photons in the event were checked in a search for invariant 
masses mi and rri2, satisfying the condition 



rrii — m n) 2 + { m 2 — m ir) 2 < 25 MeV/c 2 . In the energy region of this experiment the 
invariant mass of the pion pair in the process (Q) is less than 700 MeV. In the events, 
selected under this condition the clear u;(782) peak is seen in distribution (Fig. [18]), 
proving the domination of the process (^) in this kinematic region. The m nl parameter 
was defined as an invariant mass of the recoil photon and one of tt° mesons, closest to the 
mass of u meson. The 499 events found in the invariant mass region 750 MeV < m 7r7 < 
815 MeV were assigned to an ujtt class. 189 events with m wl outside this interval and 
> 700 MeV were assigned to the 7r°7r°7 class. Subtracting calculated contribution 
of the process ([!]) and using estimated probabilities of events misidentification for the 
processes and @ the number of the events of the process @ in the urn class were 
estimated to be equal 449. The corresponding number of events of the decay (ffl) in the 
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7r°7r°7 class is 164. The background from the process (Q) was estimated using events of the 
uj7i° class, no additional knowledge of the actual production cross section of this process 
was necessary. 

Then, for the events of the 7r°7r°7 and un classes the comparison of experimental and 
simulated distributions in ip and 9 angles was made. Here ip is an angle between the recoil 
photon and the direction of pion in the two-pion center of mass reference frame, 9 is an 
angle of the recoil photon with respect to the beam. As it is known, the distribution in 
9 for 7r°7r°7 class is proportional to 1 + cos 2 9 and is uniform in cosip. The comparison 



(Fig. |i9|a,c) shows, that in this class of events pions are actually produced in scalar state. 
On the contrary, the experimental events of the urn class (Fig. |T9|b,d) also well match 
the hypothesis of the intermediate utt state with quite different ip distribution. 

The 7r°7T invariant mass distribution of the events with m n ^ outside the 750 MeV < 
t?Vy < 815 MeV interval (Fig. |20|a) shows significant excess over background at large 
m n7T . At m n7T < 600 MeV the sum of background contributions dominates. Detection 
efficiency (Fig. 20b) for the process (|l[) was determined using MC simulation of the process 



5*7 — > 7r°7r°7 with a scalar states S, with different masses ranging from 300 to 
1000 MeV and zero width. In addition, this simulation provided information on 7r°7r° 
invariant mass resolution and event misidentification probability as a function of m^. 
After background subtraction and correction for detection efficiency the mass spectrum 



was obtained (Fig. 22). For masses in the 600 -j- 850 MeV interval the invariant mass 
resolution is equal to 12 MeV, so the 20 MeV bin size was chosen. At higher masses the 
resolution improves, reaching 7.5 MeV at 950 MeV, thus the bin size of 10 MeV was 
used for higher masses. 

For the spectrum normalization the events of the process (||) with 7 photons in the final 
state and reconstructed 3n° mesons were analyzed. All other selection criteria were the 
same as in the present analysis. The number of observed events of the process (|5]) together 



with the PDG Table value for £>(</> — > 777) = (1.26 ± 0.06)% |I I] provided independent 
estimation of the total number of <fi mesons produced. Such an approach minimizes 
systematic errors corresponding to inaccurate simulation of tails of distributions over 
parameters, used for event selection. The measured branching ratio of the decay ([!]) for 
the restricted mass range m nn > 700 MeV is equal 

B((f) -> 7r°7r°7) = (1.00 ± 0.07 ± 0.12) • 10~ 4 , (40) 

and for > 900 MeV 

B((f) -> 7r°7r°7) = (0.50 ± 0.06 ± 0.06) • 10~ 4 . (41) 

Here the first error is statistical and the second is systematic, which was estimated to be 
close to 12%. The systematic error is mainly determined by the following contributions: 



• the background subtraction error, which decreases almost linearly with the increase 
of invariant mass m n7T and is 5% on average; 

• error in the detection efficiency estimation, which increases with m nn and is equal 

to 8% on average; 

• systematic error in B((p — > 777) is equal to 5%. 
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The invariant mass spectrum was fitted with a smooth curve (Fig. |22[) according 
to Refs.|2l], and further used for simulation of the decay ([[]). As a result the detec- 
tion efficiency of the process ([!]) was estimated ~ 15% for invariant masses within the 
600 4- 1000 MeV . This efficiency was used in the fitting of the ^-resonance excitation 
curve. The cross section was described as a sum of the processes ([!]), (^), and (Q) with 
radiative corrections taken into account (Fig. ^l|). The background due to the process (§) 
was estimated by fitting the detection cross section of the events of the ujtt type. The 
background from the process (|5|) was obtained from simulation. The beam energy spread 
and the accuracy of the beam energy determination were also taken into account during 
fitting [0, n|. As a result, the following value was obtained: 



B(<p -> 7r°7r°7) = (1.14 ± 0.10 ± 0.12) • 10" 



(42) 



which, in contrast with fl40|) and (PI) , is valid for the whole mass spectrum. In the 
systematic error estimation the following considerations were taken into account. In com- 
parison with the results flip]) and (|4T|) the accuracy of normalization (3%) and efficiency 
estimation (5%) are higher here, the background subtraction error (5%) is the same, but 
an additional systematic error of (6%) exists, due to uncertainty in extrapolation of the 
invariant mass spectrum into the region m n7T < 600 MeV. Smaller systematic uncertainty 
has a ratio of branching ratios: 



7T 



7T 7 ) 



(0.90 ±0.08 ±0.07) ■ 10" 



B((f) Tjj) 

Assuming that the process ([!]) is fully determined by /o7 mechanism, using the relation 
B(f -> 7T+7T-) = 2B(/ -> 7T°7r ), and neglecting the decay -> KK'y [gl 
estimate from fl42|) 

/ (980)7) = (3.42 ± 0.30 ± 0.36) • 10" 4 . 



(43) 



we can 



Bt 



In the 7r°7r° invariant mass spectrum in Fig. 21 the /o- meson peak is clearly seen. 
The visible peak position is close to 960 MeV, with the width about 100 MeV. The 
analysis of our preliminary data together with other already known properties of the /o 



meson, performed by N.N.Achasov [p2 |, permits to assume, that the 4-quark component 
constitutes a large part of the f meson. Actually, the observed decay probability about 
10 -4 together with the invariant mass spectrum (Fig. |2l|) can be only due to s quarks 
constituting noticeable part of the /q meson. The results of the fitting of the mass 
spectrum, performed using formulas from Refs.[[2~T|, RTF] , are following: 



m f = (971 ± 6 ± 5) MeV, T f (m f ) = (188 



ttl) MeV, 



9fK+K- 
47T 



2.10 



-0.56 (jCV i 



9%< 



47T 

= 4.1 ±0.9. 



0.5li^ GeV> 



2 I 2 

9fK+K-l 9fn+ir- 

The systematic error in mass is determined by the detector resolution. Although the 
values of the constants are strongly model dependent, their ratio is almost model inde- 
pendent. The value of the coupling constant g^ KK jAm obtained from the fit agrees with 
the predictions of 4-quark model (2.3 GeV 2 PT|), but contradicts the value predicted by 

and by almost 3 standard deviations higher than molecular 
11)- ' 



2-quark model (0.3 GeV 2 |T 



model prediction: (0.6 GeV 
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Of course the more complex mechanism of the decay under study cannot be excluded, 
e.g. contribution of the heavy and broad a state fl3l| . But such a state can probably 
produce a smooth invariant mass spectrum, not masking resonance signal from the fo 
meson in the mass region m n7T > 900 MeV. To take into account all mechanisms of the 
transition one have to perform a simultaneous analysis of all fo data within the framework 
of one model. 

In conclusion we would like to emphasize that the 

7r°7r°7 radiative decay was 

observed for the first time and its branching ratio was measured. It was shown, that the 
fo (980)7 transition mechanism dominates in this decay. Invariant mass spectrum of 7r°7r° 
system and the rate of the decay give grounds to the assumption of 4-quark structure of 
the fo meson. 



6.7 Search for the — ► rjir j decay 

The — > 777r°7 decay can be described as a radiative electric dipole transition from 
a light vector meson state into a scalar state like ao(980) with emission of a photon. 
Although electric dipole transitions are common in radiative decays of heavy quarkonia, 
only few were observed in the light quark mesons and none of 0(lO2O)-meson. Since even 
lightest known scalar resonances are quite heavy, the energy yield in the decay is low 
(< 100 MeV) leading to a small phase space for a photon. A number of estimations 
exist on the — > f]TT 0/ ~f branching ratio in different models plfl. Since the values strongly 
depend on the quark structure of the scalar intermediate state, the decays of this type 
could be a unique probe of the structure of light scalar mesons. It was suggested, that the 
decay proceeds mainly through the a (980) intermediate state and the decay probability 
is determined by its structure. At present the structure of ao(980) meson it is not well 
established and several theoretical models exist, including simple two-quark model, K°K° 
molecular model, different 4- quark schemes. The values of ?77r°7 branching ratio 
calculated in these models vary in the range of 10~ 6 -v- 10~ 4 , making the — > i]7T 0/ -f decay 
a good probe of the inner structure of the a (980) state. 

Because of great significance of the — > ^7r°7 decay, two independent analyses were 
made. The results obtained are described below. 



6.7.1 Analysis 1 

The main background sources for the process under study (||D are the following 0-meson 
decays: e + e~ — > — > 7r°7r°7 — > 57 (0), e + e~ — > — > 777 — > 37r°7 -^77 (Q), e + e~ — > — > 
K$Kl — > neutrals (fD and non- resonant process e + e _ — > ujtt — > 7r°7r°7 (||). 

The expected number of events of the process (0) at a branching ratio of Br {(f) —* 
^7r°7) = 10~ 4 is about 300, while the number of background events (||) is 3- 10 4 . Although 
this process does not produce 5-7 final states, the topology of the process (0) could be 
faked due to merging of close photons and loss of soft photons through openings in the 
calorimeter. Total statistics of the process (0) is 8 ■ 10 5 K$Kl events with K$ — > 7r°7r° 
decays. The Kl-s produced in 0-meson decays are slow and their decay length is about 3 
m, while the nuclear interaction length in Nal(Tl) is about 30 cm. Characteristic feature 
of fake 5-7 events produced by K S K L decays due to either nuclear interaction of K L -s or 
their decays in flight is an energy-momentum imbalance and poor quality of at least one 
photon in the event. 
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Primary event selection was based on simple criteria ( p0|) for 57- final state. Such 
criteria greatly reduce background from the processes (|) and @, not affecting the decays 
©0), and®. 

Next step in the event selection was based on photons quality parameter ( and kine- 
matic fit parameter Xe- The requirements were imposed that ( < 0, for all photons. Then 
kinematic fit was performed under the assumption, that selected events are e + e~ — > 57 
ones and corresponding parameter Xei describing the likelihood of this assumption was 
calculated. Events with x% > 10 were also rejected. Study of simulated events of true 
57 processes (|l|, ||) and (P) shows, that the x\ an d C cu ts reject less than 15% of true 
57 events while suppressing the process (||) by a factor of 3 and making the expected 
background due to the process (^) very small, of the order of 10 events. It should be 
noted, that in contrast with other background processes the simulation of the process 
(D) is much less accurate due to nuclear interaction of K^, and data analysis may not 
exclusively rely on this estimation. 

Characteristic feature of the process (§) is that its events must contain two photon 
pairs with invariant masses of rj and 7r° mesons. Simulation shows that if the energies 
of photons in an event are enumerated in descending order in energy, the photons from 
7] — > 77 decay have the numbers of either 1 and 2 or 1 and 3. Corresponding experimental 
and simulated distributions in m 12 and m 13 with an additional requirement, that the rest 
of photons contain a pair with |my — m n o\ < 20MeV/c 2 and more stringent photon 



quality requirements ( < —5. are shown in Fig. ^3[ Background estimations are based 
on PDG table value for cf) — ► 777 decay branching ratio ||11| , our study of e + e~ — > um° 
reaction and <fi — > 7r 7r°7 decay J7[. Distribution of experimental events in mi2®m 13 shows 
an enhancement centered at ^-meson mass, while sum of background processes is nearly 
flat in this region. The sum of all simulated background processes, where each one was 
normalized to the number of events expected for a given integrated luminosity and total 
number 0-meson decays, describes the spectrum outside the enhancement quite well. If 
this enhancement is due to the decay (|2]), its branching ratio should be of the order of 
7 • 10~ 5 , but it is not possible to extract more accurate result and error estimates from 
these inclusive spectra, because of poor signal to background ratio and decay dynamics 
remaining obscure. 

Detailed study of the process (0) requires considerable suppression of background. It 
was done using kinematic fit with intermediate 7] and 7r° mesons taken into account. For 
each event rjTt ^ and 7r°7r°7 hypotheses were tried and corresponding xL-o 7 j and X^o^Cy 
calculated. To suppress the processes (p]) and @ the following requirements were imposed: 
X^irO-y < 20, X^o 7r o 7 > 20, and for additional suppression the process (0): C < — 4. With 
these requirements, the contribution from the processes (|l|) and (||), which are themselves 
relatively rare becomes negligible. 

The resulting spectra of the E imax /E , energy of the most energetic photon, are shown 
in Fig. ^4]. Since the recoil photon in the process (|5|) has a narrow spectrum close to 360 
MeV, E imax /Eo must be more than 0.7 in this decay. This can be seeing in Fig. ^4. 
While the spectrum in Fig. 24b is well reproduced by simulation of the process @ alone, 
the spectrum in Fig. |24] a shows excess of events over the simulation of (|5]). Absence 



of corresponding excess in the Fig. p4| b indicates, that additional events are true 5-7 
ones, because the cut < ( < 10 suppresses the number of such events by an order 
of magnitude. In the event sample with ( < —4, the background (|5|) still dominates at 
E-y m ax/E > 0.7, so further analysis was conducted with events where E imax /E was less 
than 0.7, where contribution of (Bl) is small. From Fig. Bla one can see, that this roughly 
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halves the detection efficiency for the process (0) and enhances its dependence on the r)7r° 
invariant mass. The detection efficiency varies from 1% at m vn o = 970 MeV/c 2 to 5% at 



2 



m V7T o = 700 MeV/c 

The distributions of the events over photon quality parameter ( are shown in Fig. 
Selection criteria here were the same except the less stringent requirement on photon 
quality: £ < 10. It can be seen, that simulation well describes the background from the 



process ([5|), while the excess of events in Fig. |25| a at low ( is compatible with existence 
of the decay (Q) with BR of the order of 10 -4 . Additional information could be obtained 
from the distribution in xLrO-v (Fig- HD- The distributions of good events with ( < —4 



and background with < ( < 10 do not contradict the simulation. The enhancement 
at low xL-o 7 f° r good events is clearly seen. The number of selected events is 25, from 
which 5 were estimated to be a background. The corresponding branching ratio was 
calculated using the f]ir invariant mass distribution of the selected events and detection 
efficiencies, obtained from simulation. The resulting value is (8.3 ±2.3) ■ 10~ 5 , where error 
is a statistical one. For such stringent selection criteria the systematic error in detection 
efficiency could be as high as 10 -4- 15%, but at present level of experimental statistics the 
statistical error dominates. 

6.7.2 Analysis 2 

In the decay — > rjTc '-/ we consider 5-7 final state with the following selection cuts : 

• A 7 = 5; N cp = 0; 0.8 < E tot /2E < 1.1; P/2E tot < 0.15; 6 min > 27°; 

• C < 0; xl < 25. 

The main background processes (|5|), (|^), ([5]) and (|I|) have 5 photons final state with two 7r°. 
For supressing this background besides the cut ( < two additional cuts were applied: 
1 - xio n o^ > 50 - event is not consistent with 7r°7r°7 hypothesis, 2 - xt n o > 30 - the 
kinematics of the event is not consistent with the urn hypothesis. All experimental events, 
that passed the selection cuts above, contained at least one 7r°-meson candidate. Among 
other three photons not included into 7r°, the pair of photons with invariant mass closest 
to that of 1] meson, was chosen (Fig. |28D. The visible excess of events close to //-meson 
mass could be considered as an evidence of the decay <fi — > t]ti j. For the background 
estimation we used events in two intervals of m 77 : 500 MeV < m 77 < 590 MeV and 
outside this interval: 455 MeV < m 77 < 500 MeV or 590 MeV < m 77 < 635 MeV. 

104 experimental events were found outside 77-mass interval and 96 simulated events 
from background processes (|^), (^). Their ratio is equal to 1.08 ± 0.12. Taking into 
account this value and the detection efficiency of 7.2%, we obtained B(cj) — > r/n ^ = 
(0.87 ±0.30) • 10" 4 . 

Then the following cuts were applied: x^o 7r o 7 > 50, x^o 7 < 7, ( < -4. Fig.|9| and 
30| show the distribution over xl^o- for the cut ( < —4, and over the parameter ( for 



the cut x^-kOj < 7 respectively. It is seen that at the small values of ( and x^-o-y there 
is a considerable excess of experimental events (78) over the expected background ( 34 
for the background process (|5|) and 4 for (|9|), while in the region ( > —4 and x^o^ > 7 
the experimental events number is in agreement with the expected background. Taking 
into account the selection efficiency of the process under study in the region ( < —4 and 
X^o 7 < 7 (4.5%) we have: 5(0 -> 7^7) = (1.1 ± 0.30) • 10" 4 . 

For more accurate background estimation, the events, satisfying the cut ( < 2 and 
xL-o 7 < 15, were divided into 4 classes (table |3|). 

The cm. energy dependence of the cross section in each of 4 classes, was approximated 
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Table 3: 





C<-4 


-4 < C < 2 


v 2 < 7 


78 


38 


Y 2 > 7 


58 


27 



by the sum of the cross section of the process under study, background process (|^) with 
coefficient K m and nonresonant background. The expected distributions over all main 
parameters ((, xL-o 7 ) i n eac h class of events were obtained from simulation. Only ( 
distribution for the process under study was derived from experimental events of the 
process (|l|). The fit results obtained are the following: 



K m = 1.26 ±0.18, 

B(<p -»■ 7?7r°7) = (0.96 ± 0.32) ■ 10~ 4 , 



(44) 



where K m is the ratio of found cross section of the process (||) to the expected one. 
Nonresonant background was found to be negligible. Fig.[H] shows the distribution over 
i]7T° invariant mass , corrected for detection efficiency dependence on the rjir invariant 
mass. The errors indicated include statistical errors and error in background estimation. 



The fitted curve |2T| is indicated by solid line. The accuracy of our measurements is not 
enough to obtain the coupling constant g aorpr . Its value was set to g aoV TT = 0.85g aoX ^ 
HJ). As a result of the fit the following parameters were obtained: 
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986±?B MeV, g 2 anKR /^ = (1.5 ± 0.5) GeV 2 . (45) 



6.8 Evidence of the — > utt° —> 7r + 7r~7r°7r° decay 

Recently in experiments with SND and CMD-2 detectors at VEPP-2M the study of <p 
meson rare decays with branching ratios of 10~ 4 -j- 10~ 5 became possible [[H], [7], [I7|| . One 
of the decays this kind of decays is OZI and G-parity forbidden <fi — > um° decay. The 
expected branching fraction of this decay is of order of 5 • 10~ 5 []33| , [3~4f . Predictions vary 



in wide limits depending on the nature of 0, u, p mixing and existence of direct <p — * 
transition. Because of the significant cross-section of the non-resonant e + e~ — > ujtt 
reaction in the vicinity of ^-resonance, the decay <p — > ujtt reveals itself as a narrow 
interference pattern in the cross section energy dependence. This allows to determine 
from the data both the real and imaginary parts of the decay amplitude. The <p — > ujtt 
decay was not observed by now. In our preliminary study [[U| only an upper limit 
of the decay probability was established at 5 ■ 10 -5 . The present work is based on full 
statistics of PHI-96 run and part of the data from PHI-97 run at the energy points 980, 
1040 and 1060 MeV. The total integrated luminosity is AL = 3.7 pb~ x , corresponding to 
6.3 • 10 6 0-mesons produced. 

Events selection. For search of the decay <fi — > ujtt we studied the cross section 
of the process e + e~ — > ujtt — > tt + tt~tt tt (|25| ) in the vicinity of the resonance. The 
final state consists of 2 charged particles and 4 photons. Unfortunately, in addition to 4 
produced photons, we found in selected events with the probability of ~20% one or two 
faked photons. So, for analysis we selected events with 2 charged particles and 4 or more 
photons. The reconstructed production point of charged particles should be within ±0.5 
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cm from the beam and ±7.5 cm from the center of the detector in the beam direction 
(the interaction region length a z is about 2 cm). 

Because of the large probability of faked photons production virtually all main <p- 
meson decays are a potential source of background in the search for <p — > ujtt decay. 
To suppress the background from the decays e + e~ — > (f> — * K + K~ (p8|), e + e~ — > — > 
KsKl, Ks — > n + 7i~ (^9|) the following selection criteria were applied: 

• a n7T < 140°; 

• Ionization losses of charged particles are close to those of minimum ionizing particles. 

The first condition suppressed events from the process (p9|), where minimum angle 
between pions from the Ks decay is 150°, and considerable number of events from the 
process (^), where charged kaons produce two collinear tracks. Part of events of the 
process (^8|) survives this cut due to decay or nuclear interaction of one of the kaons in 
the material of the beam pipe or drift chamber shell. Remaining charged kaon due to its 
low velocity (3 ~ 0.25 has high dE/dx in the SND drift chamber, allowing to reject these 
events using ir/K separation parameter ln(P K /P n ) (Fig. j32|). 

In the kinematic fit the energy-momentum balance and masses of intermediate particles 
(tt°,i],uj) were included. Three following hypotheses of the event origin were considered 
for each event: 

• The event originates from the process e + e~ — > 7r + 7r~7r°. The value of likelihood function 

IS X37T- 

• The event is from the process e + e~ — > 7r + 7r~7r°7. The photon recoiled mass m rec , was 
calculated. 

• The event is due to the process e + e~ — > 7r + 7r~7r°7r° The value of likelihood function is 
Xa-k- The recoil mass of ir° mesons was calculated and one m 3lT closest to uj meson mass 
was chosen. 

In case if number of found photons was greater than that in considered hypothesis, 
extra photons were rejected as fake photons. To do that, all combinations were considered 
and one was left with a minimum value of likelihood function. The same approach was 
applied in search for the best it candidates among all possible photon pairs in the event. 
The distribution of experimental and simulated events from the process (^5|) over the 
parameter x 3?r , is shown in Fig. [33|. One can see a considerable contribution from process 
(PT|) as a peak at low values of Fig- shows the experimental spectrum over m rec , 
where clear peak of rj meson from the reaction (|32|) is seen. To suppress the background 
from processes ( |31"D and (|32| the following cuts were used: 



• X3n > 25, 

• m rec . > 620 MeV. 

Data analysis. Fig. ^ shows the distribution of experimental and simulated events 
from the process (^5|) over x 47r . The considerable difference between the "tails" of mea- 
sured and simulated spectra is a sign of background, which survived the cuts. In Fig.p5|the 
distribution over ir° recoil mass 7713,,- is presented for experimental events with XAn < 20, 
simulated events of the process (^) without utt intermediate state (the model pirn was 
used) and simulation of the process ( p5|) with 7% addition of the process (|33p. The last 
distribution is in a good agreement with experimental data. 

For further analysis all events were divided into 4 classes: 



1. X4* < 20, \m 37T - 782| < 30; 

2. X 4n < 20, 30 < \m 3n - 782| < 60; 
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3. 20 < Xin < 40, \m 3T - 782| < 30; 

4. 20 < x*n < 40, 30 < |m 37r - 782| < 60. 

The visible cross section for each class was fitted using the following formulae: 



= eB{u -» 37i)(a + A(E - m$)) 



i _ z m^£± 



2 



(1 + 5), 



where E = 2Eq, <Jq is a nonresonant cross section of the process ( p5[ ) at i£ = m^, A 



its slope, e - the detection efficiency of the process fl25"D with X4tt < 40 and |m 3 „ 
782| < 60MeV, Z - complex interference amplitude, m^,r^ - meson mass and width, 
= m^ — E 2 —iEY^{E) - 0-meson propagator, B(uj — > 3n) = 0.888 - the branching ratio 
of the decay u — > 37r [jTYJ , 5 - radiative correction []3"5| , cr 47r - the visible cross section of the 
process (|33"D, - the visible resonant background cross section for z-class selection, «j 
and f3i - the probabilities for events from processes (|25|) or (|33|) to be found in the i-class. 

The fitting was performed for all 4 classes simultaneously. The class 1 with a small 
resonant background, was the most important for evaluation of o"o, A, Z parameters. 
Classes 2-4 were used to determine the background from 0-decays a^, 0^3, 0^4 . It was 
suggested in the fit, that for resonant background the distribution over m 37r does not 
depend on x 4?r , and for class 1 the background can be obtained from the expression: 
&(l>i = °>2 • (0^3/0^4) ■ The cross section of the resonant background a^i in the resonance 
maximum was found to be = (28 ±17) pb, what was about 4% of the visible cross 
section (E5|) in class 1. 



The contribution of the process (|33|) was determined from the ratio of nonresonant 
cross section for events in the intervals ±50 MeV and ±100 MeV of the w-meson mass. 
The value of the cross section (p3|) was found to be (6.9 ± 3.8)% of the cross section (|25|) . 
With the chosen cut ±60 MeV this contribution was suppressed down to (4.2 ± 2.3)%. 



The chosen selection criteria suppress resonant background, but decrease the detection 
efficiency of the process (p5|). In addition, the tracking in the drift chambers and nuclear 
interaction of pions are simulated imprecisely. Therefore, the detection efficiency, obtained 
from the simulation, was corrected in three ways. First, we processed the data with softer 
cuts on the parameters x^n an d m 37T , second, we studied a class of events with 1 charged 
and 4 photons outside the 0-peak, and third, we compared the efficiencies under different 
cuts for both data and simulation. The total correction obtained from experimental data 
was found to be 12%. It was mainly due to inaccuracy of the simulation of xa-it parameter 
distribution and drift chamber reconstruction errors. The corrected detection efficiency 
was 16.9 ± 1.7% at E = m^. Its value is almost constant in the energy range under study. 
The systematic error was estimated as 10%. The number of events in classes 1-4 («i, Pi) 
was obtained in part from simulation or from the fit results. The parameters «j, $ of 
the energy dependence in linear approximation were obtained also from the fit with the 
constraint, that detection efficiency is energy independent. 

The number of fit parameters, describing energy dependence of the cross sections in 4 
classes of selected events, was 12. In each class the cross section was measured in 14 energy 
points. At x 2 /d-f- — 35/44 the following values of main fit parameters were obtained: 

ctq = (8.6 ±0.9) nb, 
A = (0.090 ± 0.011) nb/MeV, (46) 
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Re(Z) = 0.108 ±0.026, 
Im(Z) = -0.127 ±0.027. 

The visible cross section in class 1 and fit are shown in Fig. |3?|. Despite 4% resonant 
background, the interference pattern is clearly seen. The interference amplitude can be 
presented in the form Z = \Z\ ■ with 

\Z\ = 0.167 ±0.027, 

i> = (-50 ±9)°. 

The branching ratio of the decay — > is obtained in the following way: 

5(0 - con ) = U -^l = (5.7±?-g) ■ 10- 5 , 

&</> 



where = l2TrB(<f) —>■ e + e )/mfj ) = 4240 nb - the cross section in the maximum [|Tl[ 



To estimate systematic errors caused by possible detector instability during lengthy 
data taking runs, we processed data from 1-3 runs separately from 4-6 runs and PHI-97 
run. It was found, that all three data samples are well described by the fit obtained above 
from the whole data. 

To investigate reliability of the obtained results (PD, we changed the cuts: for the 
parameter x 47r we used cuts 50 and 100 instead of 20 and 40, for the parameter |m 37r — 782 1 
we used values 50 and 100 instead of 30 and 60. The detection efficiency increased up 
to 20%, the resonant background in the maximum was a^x = 0.206 ± 0.069 nb, that is 
20% from the visible cross section (pip. But the interference amplitude was obtained: 
Re{Z) = 0.117 ± 0.025, Im(Z) = -0.132 ± 0.028, which does not contradict to (g§) and 
thus confirms our procedure of the background subtraction. As a final result we prefer 
(|46|), because it was obtained at lower resonant background. 

Conclusion. Obtained in this work non-resonant cross section of the process e + e _ — > 
U7T is in agreement with our old result [15| for the neutral to — > 7r°7 decay: a(e + e~ — > 



CJ7T ) = (8.7 ± 1.0 ± 0.7) nb. Our preliminary result in is by 16% lower, but this 
deviation is explained by systematic uncertainty in simulation, discussed above. The 
measured non- resonant cross section o"o exceeds in two times the expected value, where 
only p(770) — * un transition is taken into account. The agreement should be significantly 
improved with radial excitations, included into calculations. The measured interference 



amplitude Z (46) is near lower edge of theoretical predictions [34]. But in [34] the known 
radial excitations of p also were not considered. Another important remark is small value 
of measured real part of interference amplitude Re(Z), which could hardly be explained 
by well known — uj mixing model [j34| |. For example, the predicted in ^3] branching ratio 



of the decay — > lott , is 1.5 times higher than measured in this work. The interference 
amplitude 

\Z\ = 0.166 ±0.027, 

measured in this work, is of six standard deviation significance. Thus, we claim the 
existence of the decay — > um° with the branching ratio of 

5(0 -> U7i°) = (5.7^-0) . 10" 5 . 
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6.9 Study of the process 



e^e 



UJ7T 



7T 7T 7 



During the study of the reaction e + e 



7T + 7T 



(ID, the 



UJ7V 



decay 



was observed for the first time with the branching ratio 5 ■ 10 -5 |6.8| . The decay reveals 
itself as an interference wave in non-resonant cross section of the process e + e _ — > um°. In 
principle, the similar picture should be observed in neutral channel e + e~ — ► um° 
(0). Really the whole situation here looks more complicated because of other 
neutral decays like 



7r°7r°7 



A , 



/o7, e l [21 



meson 

which can create the same final state and 



interfere with the (|9]) process. The interference amplitude with the — > //V° decay is 
about 10%, which is close to the value 17%, obtained in [0. In our preceding study of the 
reaction (|25| ) fi~6| , [J], we did not observe the interference because of small statistics and 
non-resonant background. 

In the present work we studied the cross section of the process (|^) in the vicinity of <fi 
meson. Besides this process, the same 5 photon final state have processes e + e~ — > <fi — > 
7r Q 7r°7 — ► 57 (|Tp, e + e~ — > — ► K S K L — * neutrals (^), e + e~ — ► — > 77 — > 37r°7 — > 77 

(§, 

Some contribution to the background give cosmic rays and particles from collider 
beams. The events from two photon annihilation with splitting showers also mimic 5 
photon configuration, because of its large cross section. The cut imposed on the total 
energy deposition and total momentum (PU[), strongly suppressed the process (^|), the 
cosmic and beam background and considerably reduced the contribution of the process 
(EJ). Two photon events were suppressed by the cut on the the energy of two most energetic 
photons E 2 < 0.8 ■ E . 

In the kinematic fit the requirement was imposed, that two tt° are found in 5 photon 
event. The corresponding parameter for experimental events and simulation of the 



processes @ and (§) is shown at the Fig. The cut X^o^ < 40 was imposed for further 
analysis. Then the 7r° recoil mass was chosen m w , closest to to mass (Fig. |39p. The cut for 
this parameter — 782 1 < 60 was used later. The additional suppression of processes 
(01) and (BD was obtained with the parameter (, describing the transverse profile of the 



electromagnetic shower (Fig. 40). For events selection the cut £ < 25 was applied. 

Analysis of data. The events selected by the criteria, described above, were divided 
into 7 following classes: 



X^o„-o 7 < 20, 
X^o„-o 7 < 20, 
20 < xk« 7 < 4 0, 
20 < xlo^ < 40, 
Xlo^ < 40, 
X^o„.o 7 < 20, 
X^o„-o 7 < 20, 



\m u - 782| < 30, 
30 < \m w - 782| < 60, 
\m u - 782| < 30, 
30 < \m u - 782| < 60, 
|m w - 782| < 60, 
7821 < 30, 



7821 < 30, 



C<-5 
C<-5 
C<-5 
C<-5; 
-5 < C < 5; 
5 < C < 25; 
AL = 6. 



The 6 photon events were put into the last class. It was done to investigate, for example, 
a possible background from beams, where sixth photon is superimposed on the event. 
In the fit such a photon was supposed to be spare. In the Table f| there are shown 
the probabilities for the process under study to be found in classes, described above. 
The probabilities for events of the process under study to be found in classes 1-7 are 
shown in Table § as well as the background resonant cross sections, normalized by the 
e + e~ — » (f) — > cross section. The process ([!]) was simulated for / 7 intermediate 
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Table 4: The probabilities oij for events of (^j) to be found in z-class, the ratio fa of non- 
resonant cross section in i-class to the e + e~ — > <fi — > rj'y process cross section, obtained by 
simulation and from the cross section fit. 



ass 


Cv- PTTi 


a; MC 


Bi MC(lTY fr>ry KsKr) • 10 3 


R- PTTI • 1 3 


1 


V 7 . 'J i l 1 — 1_ \ ' . \f 1 


0.633 


18+0 20+0 00—0 38 


0.25 


2 


0.061+ 0.013 


0.067 


0.20+0.18+0.00=0.38 


0.23+ 0.11 


3 


0.069± 0.015 


0.049 


0.43+0.02+0.00=0.45 


0.49+ 0.12 


4 


0.020+ 0.009 


0.012 


0.37+0.02+0.00=0.39 


0.36+ 0.08 


5 


0.196± 0.023 


0.202 


1.33+0.11+0.09=1.53 


1.74+ 0.23 


6 


0.042± 0.012 


0.033 


0.30+0.01+0.00=0.31 


0.40+ 0.10 


7 


0.021± 0.012 


0.006 




1.87+ 0.16 



state with the branching fraction of 0.9 • 10~ 4 . The simulation shows that the signal to 
noise ratio (S/N) in class 1 is the largest, which was used later for investigation of the 
process (|9]). The contribution of the process (|9]) in the class 2 is about 10 times lower, 
which allows to extract from the data the resonant background from the processes ([!]), 
®> ©• Using ratio of the resonant background in classes 6 and 7 obtained from the 
simulation, one can estimate the background in class 1. The ratio of the background 
processes in the regions \m u — 782 1 < 30 and 30 < \m w — 782 1 < 60 weakly depends on 
the cuts in X^-o 7r o 7 and ( parameters. For the process ([TJ) this ratio varies from 1.09 to 
1.12, while for the process @ — from 0.9 to 1.07. But in the latter case this range is 
determined by low simulation statistics of the process (|5|). The coefficient 1.1 ± 0.2 was 
taken for background recalculation from class 2 to class 1. The classes 3-7 were used 
for estimation of systematic errors, depending on different selection cuts. 
The visible cross section in each class was presented in the following form: 



i - z m > 



2 (l + £), 



where do - non- resonant cross section of the process e + e _ — > ujtt — > 7r°7r°7 at E — m^, 
A - the slope of the cross section, e - detection efficiency of the process (f|), Z - complex 
interference amplitude, m<j>, - mass and width of cf> meson, = mi — E 2 — iET^(E) 
- 0-meson propagator, B(uj — > 7r°7) = 0.085 - the decay uj — ► 7r°7 branching ratio fl~l"H , 
5 - radiative correction p5| , cr^-^ - the cross section of the process e + e~ — > (f> — > 777. 
The parameters oij mean the probability for the process @ to be found in the class i, 
the parameters fa - non-resonant cross section in the class i, normalized to the process 
e + e~ — > — > 777 cross section. The <fi meson excitation curve was described by the main 
background process e + e~ — > <p — > 777. The differences in excitation curves for different 
processes (^), (Q),([l]) were neglected at the present level of accuracy. The fitting of the 
visible cross section was done for all four classes. The parameters <7o, A, Re(Z), Im(Z), aj, 
fa except fa «i were set free, fa was found from the expression fa = (l.l±0.2)-/3 2 , «i from 
normalization ^Z^i = 1- The detection efficiency e=39%, obtained from the simulation, 
does not depend on the energy. The coefficients q.{ are also energy independent. The total 
number of fit parameters is 16. In each class the cross section was measured in 15 points. 
The following values of fit parameters were obtained (x 2 /d.f. = 63/89): 

do = (0.65 ± 0.04) 716, 
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A = (0.0065 ± 0.0018) nb/MeV, (47) 
Re(Z) = 0.036 ±0.052, 
Im(Z) = -0.186 ±0.063. 

One can see from the Table ^ rather good agreement between experimental and simulated 
parameters a iy fy. This allowed to put the total systematic error to 6% for a . For other 
parameters in ( f4"T| ) the statistical error is much higher than the systematic one. The 
visible cross section for class 1 and fitting curve with x 2 /d-f- = 1 1-9/11 are shown in 
Fig. fy]. Fitting curves for the visible cross section of the process (|9|) and for non-resonant 



background are also shown there. One could see, that in spite of imposed strong cuts in 
class 1, the resonant background is about one third of interference amplitude wave. It 
remains the dominant source of systematic error in Z. 

Conclusions. The measured non-resonant cross section of the process e + e~ — > ujtt 



(7.7 ± 0.5 ± 0.5) nb agrees with the result (8.7 ± 1.0 ± 0.7) nb from ]TJ| and with the 
result from |7j in the decay channel with charged pions e + e~ — > ujtt — > 7r + 7r~7r°7r°. The 
measured interference amplitude is three standard deviations above zero. The measured 
in |7| interference amplitude in the <fi — > ujtt decay is Re(Z) = 108 ± 0.026, Im(Z) = 
—0.127 ± 0.027. Our calculation of the interference amplitude for the decay <p — > pn — > 
7r°7r°7 gives Re(Z) = —0.079, Im(Z) = —0.053. The sum of these contributions Re(Z) = 
0.029, Im(Z) = —0.180 agrees with our measurement. We did not make a theoretical 
estimations of the contributions of the processes — > /o7> e l ~* vr°7r°7 into interference 
amplitude. Finally we list parameters of the e + e~ — > ujtt — > 7r°7r°7 cross section, obtained 
in the present work: 

(x = (0.65 ± 0.04 ± 0.04) nb, 
A = (0.0065 ± 0.0018) nb/MeV, 
Re(Z) = 0.036 ±0.052, 
Im(Z) = -0.186 ±0.063. 

6.10 Ks — > 37T° decay search 

This decay is a pure CP-violating one, similar to those of — > tttt and unlike Ks — > 
tt + tt~tt° mode, where C-conserving amplitude could be present. Considering the model, 
where CP-violation originates from the mixing of CP-even and CP-odd states K° and 
K® , it leads to physical state Ks — K® + sK®, and the branching of Ks depending on e: 
B(Ks — > 37t°) ~ \e\ 2 ■ ^B(K L — > 3tt°) ~ 10 -9 . Very small magnitude of the branching 
ratio and experimental problems with separation between this decay and background from 
the process Ks — > 27r°, which has a 9 orders of magnitude higher probability, result in 
rather high existing experimental upper limit: B(Ks — > 37r°) < 3.7 • 10 -5 . 

Analysis. In the present work the (f) — > KsKl decay with a total = 7.5 ■ 10 6 is 
used as a source of Ks- Our analysis is applied to those events, where Kl leaves the 
detector without either interaction or decay. 10 4 simulated events of Ks — > 3tt° and 
10 4 simulated events of Ks — > 2tt° were processed in the same way. Because the initial 
decay <fi — > KsKl was not simulated, the calculations take into account the branching 



ratio B{(j> K S K L ) = 0.343 ± 0.007 |TJ and the probability W Kl for K L to leave the 
detector undetected. The latter value cannot be obtained from simulation due to its low 
accuracy in nuclear interactions of Kl mesons, so it was estimated to be 18 ± 1% from 
the experimental data on — >• KsKl, Ks — > 2tt° decay. 
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The following criteria were applied to the experimental data and simulation sample 
for Ks — > 37T° decay search: 

• iV 7 > 6; N cp = 0; positive neutral trigger; 0.38 < E tot /2E < 0.5; 

• 0.15 < P/E tot < 0.3; E imin > 30 MeV; ( < 10. 

The hits in the first two layers of the calorimeter are required. The kinematic fit was done 
for the selected events, in assumption of Ks — > 37r° kinematics. The output parameters 
within the limits shown below were the following: X 2 k s ^3w° < ^ f° r ^ ne found combination 
satisfying Ks — > 37r° — > 67 decay hypothesis; minimal angle between photons of the Ks 
decay and beam axis a^ s ^ 37r o > 30°; polar angle of reconstructed K s 30 < "&K s -*3n° < 150; 
raw invariant masses of three 7r°'s: 110 < m w o 2 K s -*3ir° < 160 MeU. 

For the selection of the process Ks — > 2tt° the following criteria were applied: 

• iV 7 = 4; positive neutral trigger; 0.38 < E tot /2E < 1.1; P/£7 tot < 0.6; 

• E ymin > 20 MeV; £ 7ma:E < 300 MeV. 

The cut 0.05 < E tot /2E - P/E tot < 0.4 selects a slanting stripe in the (E tot /2E , P/E tot ) 
plane, containing most of the (ft — > KsKl events with tagged Kl (Fig. |B0|). The kinematic 
fit has been carried out also in assumption of the Ks — > 27r° event kinematics. The 
selection criteria after the kinematic fit were the following: Xk ->2n° < ^5 for the found 
combination of photons satisfying Ks — > 2tt° — ► 47 decay hypothesis; minimal angle 
between photons of the Ks decay and the beam axis ceK s ->2n° > 27.5°; the raw invariant 
masses of the two tt° mesons 110 < m^o 2> k s -*2w° < 160; the reconstructed K s momentum 
50 < Pk s ^2w" < 200; the reconstructed impact parameter —2.5 < Z Ks ^ 2 w° < 2.5. 

The results. One candidate event of the Kg — ► 3tt° decay was found. The detec- 
tion efficiency of this process is e Ks ^ 37T o ~ 14%. The detection efficiency of the process 
Ks — > 2tt° is £x s -»27T° — 46%. The total statistics of found events of this process after all 
cuts is 6869 events. 

Using the formulae 

N Ks ^o = ■ 5(0 -> K S K L ) ■ W Kl ■ B(K S -> 2tt°) • e^^, (48) 

N Ks _> 3n o = iV 03 • B((f> -> ir s K L ) ■ ^ • 5(^5 ^ 37T°) • £ M .o, (49) 

where iV^j ^3 are the numbers of produced </> mesons, we found: 

S(K 5 -> 3tt°) = 1.8 ■ 10" 5 

the corresponding upper limit was evaluated to be 

5(^5 -> 3tt°) < 6.9 ■ 10" 5 , C.L. 

Our result is close to the table limit of < 3.7 ■ 10 -5 [|TT|. 



6.11 Study of e + e — > e + e 7 and e + e — >• e + e 77 processes 



The process e + e~ — >■ e + e~7 (|26| ) is a QED process of the order of a 3 . Its study is 
interesting not only for QED testing, but also because it is important as an integrated 
luminosity check. It is also a source of background in searches for rare hadron processes. 
To compare experimental results with the QED the simulation was used from the work 



The process e + e — > e + e 77 fl2?D is a QED process of the 4-th order in a. Its 
investigation is also useful as a QED test, especially because it was studied in only one 
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experiment [|15j with ND detector, where 223 events were seen. The measurement of its 
cross-section is important for study of <fi — > ne + e~ n — ► e + e~7 decays, where this process 
is a main source of background. 

The study of the processes is based on the data sample with a total integrated lu- 
minosity of 4.3p6 _1 in the energy range 985 4- 1040 MeV, while the production cross 
section measurements are based on the data subset, corresponding to ~ lpb' 1 integrated 
luminosity. 

Analysis of e + e~ — > e + e~7 (|26|) . To study the process (|26|) the following selection 
criteria were used: 

• positive charged trigger; N cp = 2; 1 < iV 7 < 3; Ri, R 2 < 0.5cm; \Zi\, \Z 2 \ < 10cm; 
9 min > 36°; |A0 ee | > 5°; 

• no hits in muon system; E tot /2E > 0.8, P/E tot < 0.15; Xe < 15. 

Extra photons were allowed to prevent the events loss due to accidental signals in calorime- 
ter due to beam background. For events passing these selection criteria, the experimental 



spectra and QED simulation |3(J are shown in Fig. f|2[ The events number in simulation 
corresponds to 190 nb~ x integrated luminosity. The cross section values were obtained 
separately for each scan. Resulting values are lieing between 15.5 17.5 nb, giving an 
average cross section a = 16.7 ± 0.1 nb , with x 2 — 24/5d.f.. After that all scans were 
processed again with the same configuration of trigger, the event losses due to trigger inef- 
ficiency was taken into account and luminosity based on e + e~ events was used. As a result 
the difference of measured cross sections in different scans was reduced to 17.1 -f- 17.9 nb, 
that is less than 5% . The average cross section was found to be a = 17.6 ± 0.1 nb , with 
X 2 = 7/5. But the simulated cross-section was significantly higher a S i m , = 19.5 ± 0.3 nb, 
so the difference between the experiment and simulation was ~ 10%. It was assumed 
that this difference appears because of different angle 9 cuts in our search and in luminos- 
ity calculations (in the first case angle is determined using drift chamber, in the second 
- using calorimeter). To study this problem in two scans PHL9603 and PHL9604 the 
collinear events were processed. The "e + e~" events were selected with the same criteria 
on 9, E tot /2Eo, P/E to t as events "e + e~7". Then the value of N eei /N ee was calculated for 
experiment and simulation. The results are shown below: 

iV ee7 /iVe e (PHL9603) = 0.0138 ± 0.0002, 

N ee7 /A / e e (PHI_9604) = 0.0141 ± 0.0001, 
N eei /N ee (Simulation) = 0.0142 ± 0.0003. 

One can see that in this case the results in both scans are in good agreement with each 
other and with simulation. For 2 scans the integrated luminosities were recalculated and 
the visible cross section was obtained: a exp , = 19.0 ± 0.2 nb, a S i m , = 19.5 ± 0.3 nb. 
Taking into account the detection efficiency 80.5%, we obtained the cross section in the 
polar angle larger than 36° and acollinearity angle larger than 5°: o exp . = 23.6 ± 0.3 nb, 
Csim. — 24.2 ± 0.4 nb. The systematic error is determined mainly by the error in the 
integrated luminosity and is equal to 5%. The simulation and experiment are in good 
agreement and the statistical error is at a level of 1%. So, to increase the accuracy it is 
necessary to reduce the systematic error. 

Analysis of the process e + e~ — * e + e~77(|27|). To study the process (|27| ) the follow- 
ing selection criteria were used: 

• charged trigger N cp — 2; 2 < A 7 < 3; Ri, R 2 < 0.5 cm; \Z\\, \Z 2 \ < 10 cm; 9 min > 36°; 
|A</>ee|>5°; 
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• no hits in muon system; E tot /2E > 0.8, P/E tot < 0.15; 

• to cut events with tt° decays, the invariant mass of two photons m 77 was required to be 
outside the interval from 110 to 170 MeV; £ 7 min ) > 50 MeV; Xe < 15 - 



For selected events the spectra are shown in Fig.[43] together with simulation ||38||. The 
number of events in simulation corresponds to integrated luminosity of 4297 nb -1 . The 
main background for this process comes from the reaction e + e~ — > e + e~7 (|26"D with 
additional soft photon from the beam background. To suppress these events the threshold 
in minimum photon energy was placed at 50 MeV. For this process the ratio N eeil /N ( 
in separate experimental runs and simulation was obtained: 



ee 



iV ee77 /iV ee (PHI_9603) = (0.10 ± 0.01) • 10~ 3 , 

A^ ee77 /iV ee (PHI_9604) = (0.12 ± 0.01) • 10~ 3 , 

N eni /N ee (Simulation) = (0.10 ± 0.01) • 10~ 3 . 

Using corrected luminosity the visible cross section was obtained: a exp , = 0.14±0.03 nb, <J s im 
0.135 ± 0.006 nb. Taking into account the detection efficiency 32.6%, the cross section 
for polar angle greater than 36° and acollinearity angle larger than 5° was obtained: 
&exp. — 0.43 ± 0.09 nb, and cr s i m . = 0.42 ± 0.02 nb. The systematic error like in previous 
analysis is determined mainly by error in luminosity and is equal approximately to 5%. 
As one can see, the simulation and experiment are in good agreement, but the statistical 
error in experiment is at a level of 20%. So it is possible to improve accuracy with more 
statistics. The energy dependence of cross section (p6| ) and (p7|) is 1/E 2 in agreement 
with QED predictions. 

6.12 Study of QED process e + e — ► 37 

Earlier the three photon annihilation process was studied in high energy region 2E ~ 
10 GeV [^, In our previous experiment with ND detector this process was studied 
as a QED background [ 15 1 . Investigation of (|ll|) process at low energy near 0-meson peak 



has an advantage, because its cross-section energy dependence ~ a 3 /Eq provides us with 
high statistics. This gives us a possibility for more detailed QED test in the third order of 
perturbation theory. Besides, the e + e~ — > 37 process is important as a background source 
for 0-meson radiative decays e + e~ — ► </> — > 777 — ► 37 (|7]) and e + e~ — ► (p —* tt°j ^37 

In the present study, the experimental data of PHL9603 scan were analyzed. The 
total integrated luminosity is AL = 822 nb~ x . The following selection cuts were applied: 
Nry = 3, 6 min > 27°, 0.7 < E tot /2E < 1.1, x% < The main background comes from 
the processes (|7|), (g) and e + e~ — > 77 (|TDj) . Taking into account presence of a quasi- 
monochromatic recoil photons in radiative decays, we imposed the cut E ltmax /E < 0.974 
to suppress (§) process. To suppress (|7|) process, the photons with energy in the 0.65 < 
Ej/Eq < 0.77 interval were rejected. The background from (|T0| ) process, was significantly 
reduced by condition imposed on acollinearity angle in the azimuth direction between 
two most energetic photons |A0 12 | > 6°. The detection efficiency was obtained by MC 



simulation, and for described selection criteria it equals (10.3 ± 0.3)%. Fig. 44 shows the 
energy dependence of visible cross-section of the e + e~ — > 37 process. The value of the 
cross section at Eq = 510 MeV was found to be: a exp , = (1.82 ± 0.21 ± 0.25) nb. The 
first error is statistical, the second one - systematic. The simulation with the standard 
QED matrix element gives: u S im. — (2.01 ± 0.05) nb. The cross-section of background 
processes at peak does not exceed 5% of that of the process (|TT|) . The value ~ 5% is 
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due to resonant background. The detection efficiency for all background processes (|7|), 
(||), (|10"D , obtained from simulation was not higher than 0.25%. 



Figs.H5|, 06], 07] show the softest photon energy spectrum, distributions over photon 



pair invariant mass and acollinearity angle in the azimuthal direction; points correspond 
to the experimental data, histogram - to the MC simulation. These distributions show 
good agreement between experiment and simulation, confirming the validity of QED at 
our level of accuracy. 



7 Physical results from MHAD-97 experiment 



7.1 Investigation of e 



7T + 7T 7T 7r , 7T + 7T 7T + 7T 



reactions 



Processes of e + e~-annihilation into four pions attract attention due to the following rea- 
sons. First, in the 2Eq = 1 2 GeV energy region these processes dominate and deter- 
mine the main part of the hadronic contributions into the anomalous magnetic moment of 
muon and into QCD sum rules. Second, the processes e + e~ — > (p — ► 7r + 7r~7r + 7r~ (^4j) and 
e + e~ — ► 7r + 7r~7r°7r° (^) are an important source of information for hadron spectroscopy, 
in particular, for study p-meson radial excitations. 

According to the existing data, in the 2E = 1 2 GeV energy region there exist two 
radial excitations of the p meson: p(1450) and p(1700) (TT|. Determination of parameters 
of these states and their interference with the p meson, can be found, for example, in work 
||41|| , |02fl . A possible mixing of these excited p-states with the exotic ones (for example, 
4-quark states) is discussed in |4"3||,P4"|. Recently some experimental evidence appeared in 
favor of p x (1300) state existence f45fl , which possibly is not a conventional quark-antiquark 
meson |06|| . 

In the past these processes were studied at VEPP-2M 0,|47], H, DCI[[||, ADONEg] 
e + e~ colliders. Statistical accuracy achieved in these experiments is ~ 5%, with a sys- 
tematic error of ~ 15%, and the discrepancy between different experiments is sometimes 
as large as ~ 20% [15, 0_7J. Therefore, measurements with smaller systematic errors are 
needed to clarify situation with multi-hadron production. 



7.1.1 e + e — > 7r + 7r 7r + 7r process 

To select events of the process e + e~ — > 7t + 7t~7t + tt~ ([24]) the following selection criteria 
were applied: 

• N cp > 4, iV 7 > 0, N wire < 30, E tot /2E < 0.8. The energy dependence of visible cross 
section, obtained under these conditions, is shown in Fig. ^S] with error bars indicating 
only statistical errors. The following sources of systematic errors were taken into account: 
charged particles reconstruction errors, systematic errors caused by the use of the selection 
cuts, i.g. cuts on the recalculated trigger and the number of hit wires, inaccuracy in 
luminosity determination. As a result, the systematic error was estimated as ~ 15% 
(preliminary). The detection efficiency was determined from the Lorentz- invariant phase 
space simulation (LIPS). As it is seen from two-pion invariant mass distribution in Fig. 
the experimental distribution has a bump around 750 MeV, which is absent in LIPS 
simulation. This can be considered as a manifestation of p meson in the four charged 
pions final state. In future, simulation with the pirn intermediate state will be made to 
compare experimental data with the simulation. Refinement of the systematic error and 
the use of the whole MHAD9702 experiment statistics is also planned. 
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7.1.2 e + e — > 7r + 7r 7r°7r° process 



The process ([33]) was studied in the final state with 2 charged particles and 4 photons. 
In the energy range above <fi meson this process dominates. Background comes from 
e + e~ — > K + K~ process, from QED processes e + e~ — ► e + e~e + e~, e + e~77, and also some 
background comes from cosmic rays and collider beams. To select events of the reaction 
( |3~3|) , the following criteria were applied: 

• N cp = 2; iV 7 > 4; E np /2E > 0.3; R u R 2 < 1 cm; Z 1; Z 2 < 10 cm; |Z X - Z 2 | < 5 cm; 
The calculated trigger corresponds to the experimental trigger for charged particles with 
refined thresholds. 32195 events survived the cuts. 

• For energies E np < 550 MeV a background from the e + e~ — * K + K~ process is sig- 
nificant. To suppress this background, additional cuts were used, related to ionization 
losses in drift chambers and acollinearity angle: a v + n - < 160°, dE/dx n ± + dE/dx^T < 
5 ■ dE I dx m _i_p_. 

• The kinematic fit was made under assumption, that charged particles are tt mesons and 
4 photons originate from two tt°. To suppress the remaining background, the following 
cuts were applied: Xe < 50, xl+^-^^o < 50, N n o = 2. 

Data analysis. The detection efficiency for the ([33]) process, obtained from simula- 
tion, is practically energy independent and equals to e — 0.30 ±0.01 Various intermediate 
states such as lutt°, /2P , afir^ and others can contribute to the reaction under study. 
The 7r°-meson recoil mass distribution in Fig.[50] shows that it is possible to separate the 
(f25|) process. It is not possible to separate other intermediate states, because they are 
very broad and distorted by interference effects. 



To separate the ( |25|) process, the selected 4tt events were divided into two classes 

• 750 MeV < M rec ^o < 820 MeV ( cott , class 1), 

• M rec , n0 < 750 MeV or M rec _ n o > 820 MeV (not cjtt , class 2). 

Using the detection efficiencies obtained from simulation with the un intermediate 
state for the class 1, and with Lorenz Invariant Phase Space (LIPS) model for the class 
2, the total cross section of the ( p5|) process was calculated. In Fig.|51| the cross section 
of the process (|25|) is shown together with the process ([33]) with all other intermediate 
states without interference between them. Fig.|52] shows the cross section of the process 
(p5|), in comparison with existing experimental results. It is seen that the measured cross 



section is lower that obtained in ND experiment |L5[], but in good agreement with earlier 
OLYA detector data |||. 

Simulations of different intermediate states like um°, p 7r°7r° and LIPS shows, that the 
detection efficiency varies within 5% range, which might be a source of systematic error. 
Another sources of systematic error are errors in event reconstruction and inaccuracy in 
simulation of 7r-meson nuclear interaction. In total, the systematic error is estimated still 
to be about 15%. 

7.2 The study of the process e + e — ► tt + tt ~7r° 

The process e + e~ — > tt + 7t~7t° (|30|) in the energy range 2E =1.1-1A GeV is interesting 
for several reasons. It is well known, that near u and <fi resonances the cross section 
e + e~ — » u, <f) —> pit — >• 7T + 7r _ 7r° is well described by the Vector Dominance Model (VDM). 
New precise measurements in non-resonant region will allow to investigate the limitations 
of VDM and determine possible contribution from heavier states like u;(1420) or cj(1600) 
These states, decaying into 7r + 7r~7r° are considered now as radial excitations of 
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u(782). Due to existence of cu-meson decay into 7r + 7r - with a 2% probability, the process 
(j30|) can proceed via e + e~ —>■ p — > c<J7r — ► 7r + 7r~7r° mechanism. As a result, the p- 
u; interference can be observed in 7r + 7r~ mass spectrum. According to ||51|| , the radial 
excitations p(1450), p(1700) (Tl|] can increase the p-u interference significantly. 

The process fl30|) was studied earlier in the energy range up to 1.4 GeV with ND 



detector at VEPP-2M |T5[]. The results and the VDM prediction are shown in Fig. [55|. 
One can see, that the measured cross section is significantly higher than the predicted 
one. 

In analysis of the process (^) the events with two charged particles and two photons 
were selected. To suppress the cosmic and beam background and contribution from the 
process e + e~ — > K S K\ the cuts were imposed on the location of the production point 
with respect to collision center: Ri,R 2 < 0.3 cm, Z\ } Z 2 < 6 cm. The main background 



processes for (|30|) are 



e+e-^e+e-(7) (50) 

e + e~ — > 7r + 7r~(7) (51) 

e + e~ — > e + e~77 (52) 

e + e~ -> tt+tt^tt (53) 



The process ( p2|) has the same final state as (|30|) : The processes ( pOjplD can mimic the 
process under study ( pOf ) because of the showers splitting in the calorimeter. To reduce 
the contribution of the process (|50|) the cut E 1 ma x/E < 0.8 was imposed. Because the 



processes (pO|, 3l , p^) have collinear tracks in drift chamber, the cut a w + w - < 160° was 
added to suppress their contributions. The main background comes from (p3|), where 
two photons from 7r°-decays are lost. We applied the cut a 77 < 100°, to decrease the 
contribution of this process. Then we applied kinematic fit to the selected events. To 
suppress the process (^3|), we used strong cut Xe < 10- The two-photon invariant mass 
m 77 distribution for events which survived the latter cut, is shown in Fig. For further 
analysis the events with m 77 close to 7r° mass within the interval 105 MeV < m 77 < 
165 MeV were selected. The cross section can be described as a sum of two terms: 

cr e+e -_ >7r+7r - 7r o(s) = a(s) + cr 7 (s), (54) 



where a is a cross section ( |30D with photon energy less than 40 MeV, cr 7 - is the cross 
section with photon energy larger than 40 MeV. The latter value, obtained from known 
cross section of the processes e + e~ — > u, (f) — > 3n, can be calculated with the accuracy of 
5%. The value a(s) is the investigated cross section with the radiative correction not more 
than 10%. The detection efficiency of the process with E 1 > 40 MeV sharply decreases 
with energy, because the selection criteria suppress events with hard photons, emitted by 
initial particles. 

The total cross section (501) was presented in VDM form: 



1-4 



2 



a(s) = Un^L, (55) 



.4 



v / r y ^ 37r (s)r y ^ ee m^e i 
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where M, V, 9 — mass, width and relative phase of vector meson respectively. The values 
To;',w"- > 37r(s)r w / ia; //^ ee and 0u/,u>" are free fit parameters. Other values were taken from the 
Tables |jTT[ . The fit results are the following: 

B^KBuS^ee = (0.14 ± 0.02) • 10- 4 
B^ 37r B^ ee = (0.46 ± 0.23) • 10~ 7 
9^ = 10° ± 10°, 9^ = 170° ± 24°. 

Disagreement with PDG data might indicate the existence of additional contributions 
besides u, <f), u/, 10" . In particular, in the cross section Q55D the interference with processes 
e + e _ — > p,p' — > urn — > 37r was not taken into account. Nevertheless, the function, 
which was used, describes the measured cross section rather well. This allows to account 
for radiative corrections and to calculate the cross section in Born approximation. The 
values of Born cross section in different energy points were calculated in the following 
way: o = ((7™p. — &™ s ')/5e. The results are shown in table [|. 



Table 5: Experimental cross section a e + e -^+^-^o. 



Total energy 


Cross section 


Total energy 


Cross section 


2E , MeV 




2E , MeV 


cr e + e -^ 7r + 7r - 7r o, nb 


1080 


2.41 ± 0.44 


1220 


4.36 ±0.41 


1100 


2.71 ±0.31 


1240 


3.56 ±0.35 


1120 


2.78 ± 0.49 


1260 


3.62 ±0.34 


1140 


2.69 ±0.49 


1280 


3.66 ±0.28 


1160 


3.66 ±0.43 


1300 


2.89 ±0.25 


1180 


3.22 ±0.40 


1340 


2.99 ±0.26 


1200 


4.05 ± 0.42 


1380 


3.04 ±0.18 



The systematic error in the cross section is about 12%. It is determined by the 
detection efficiency error (10%), background subtraction error (5%), error in luminosity 
(5%). Fig. [57] shows good agreement between old ND data || and new more accurate 
results, obtained in the present work. 

To study intermediate state in (j3(j), we measured invariant masses of 7r-meson pairs 
in the final 37r state. The intermediate state might be pn and much less probable 
with the decay ui — > 2tt. The interference between these two intermediate states can be 
observed if one compare mass spectra of tt + tt~ and n ^. Experimental data in Fig. [58 



show clear peak in ir + n~ mass spectrum, which proves the existence of p-u interference 
in 37r final state. The relative number of events in the peak area in Fig. 58 leads to the 
result 

N(U -> TT + TT") /N(p -> 7T + 7T-) = 0.047 ± 0.011. 

which is 4 standard deviations above zero. To obtain the p — u interference phase, the 
simulation was done with different values of the phase. The experimental data agrees the 
best with the phase value of 0°, which is predicted in VDM model. 
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7.3 Study of the reaction e + e — > K$K L 



The cross-section of the process e + e — ► K$Kl (|T3|) is known with a high accuracy only 
in the energy range close to the 0-meson peak. In 1982 the reaction (|TBp was studied with 



the DM1 detector (Orsay) in the energy range 2E = 1400 4-2200 MeV ]52|. At the same 
time this process was measured in the energy range 2E = 1060 4- 1400 MeV with OLYA 
detector (Novosibirsk) |53| . In both experiments the achieved accuracy was not high, so 
new measurements are desirable. 

In the present work we analyzed a part of statistics of the scan MHAD_9702 with 
integrated luminosity of ~ 1.8 pb^ 1 . The process e + e~ — > K S K L was studied in the decay 
mode Ks — > 7r°7r° (|T7|). The main background comes from the non-resonant process 
e + e~ — > u)-k q — > 7r°7r°7 (Q). Cosmic and beam background are also present. To select 
events of the reaction (|T7]) the following cuts were applied: 

• N y > 4, N cp = 0. 

• Cosmic events were rejected using SND calorimeter for reconstruction of cosmic muon 
tracks. The procedure was tested on simulated events of the process ( |i~3"|) at 2E = 
1100 MeV. From 1733 simulated events only 3 events were erraneously recognized as 
cosmic tracks, while in the experimental data 1716 such tracks were found among 2319 
events. 

• The remaining events should have two 2tt°. The invariant masses of photon pairs must 
be within the range: 115 < m gammai < 155 MeV. 

• Invariant mass of 7r°-meson pair 400 < m^o^o < 600 MeV is compatible with Ks mass. 
The cut ( < was used to suppress possible background from the Kl decay or its 
nuclear interaction in the detector material. Events passing these cuts belong either to the 
process fll3"|) or to the process e + e~ — > KsK^ with emission of photon by initial particles 
and subsequent return back to the (f) peak (|14"D, (§J). At 2E > the radiative photon in 
the process (|T^) carries the energy up to ~ 300 MeV, leading to a significant difference 
in kinematics with respect to the process (|T3|) . In particular, the energy of Ks meson 
decreases and as a result the detection efficiency with the additional cut Ek s /E > 0.94 
goes down. The cut Ek s /E < 1.06 rejects significantly the process (H). The rest of the 
events of the process (|J) is rejected by cuts E tot /2E > 0.8 and P/E tot < 0.15. Beam 
background events have photons located mainly in the small polar angle region. So, the 
cut 6 1 m in(K s ) > 31.5° was applied. 108 experimental events passed all cuts listed above. 
Distribution over m^o, invariant mass of 7r°-mesons pairs is shown in Fig. |53|. Clear peak 
at Ks mass is seen. The detection efficiency to the investigated process (0) depends on 
E Q , decreasing from 3.4% at 2E = 1100 MeV to 1.8% at 2E = 1380 MeV. Efficiency 
to the process ( |i~4"D varies from 1.9% to 0.4% in the same energy interval. Efficiency of 
the process (01) is < 0.1%. 



The visible cross section of the process ([13| ) is shown in Fig. [54] as well as expected 
cross section, calculated taking into account 0-meson resonance and radiative corrections. 
Experimental data from the OLYA detector are also shown in Fig. |54j. One can see, that 
the data obtained by both detectors are in good agreement and within experimental errors 
do not contradict theoretical estimations. 

7.4 Search for the process e + e — > KsKlk® 

KKn state could appear in e + e~ — > 07r° and e + e _ — > KK* processes. The search for 
the process e + e~ 4>tt° was carried out earlier in ND experiment in the decay mode 
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K S K L -> 2n°K L 0, H, in CMD experiment and in DM1 experiment at DCI 
| 55fl in the process e + e~ — > 07r° — > K + K~n Q . The study of the process e + e~ — > KK* 



in the channel KK* — > KsK ±r K T was conducted in DM1-DCI experiment |56| in the 
energy range 1400 ^ 2180 MeV and in the energy range 1350 2400 MeV in DM2-DCI 
experiment f57|. The channel e + e~ — > XiT* — > K + K~tt° was studied in DM2-DCI as 



well. 

New measurements, carried out with SND detector, can significantly improve our 
knowledge of these processes. In particular, new measurements could clarify the nature 
of C(1480) state [58|, found in IHEP (Protvino). With quantum numbers J PC = 1 , 



and decay into 07r°, this state is a possible candidate for exotic hybrid or 4-quark state. 
The model exists, where C(1480) is identical to p(1450)-meson f5"9fl . The electron width 



of C-state depends on its structure, so the study of C(1480) production in e + e~ collisions 
could reveal its contents. Although the C mass is larger than VEPP-2M maximum center 
of mass energy, its left slope could be observed with SND in the process e + e~ —>■ </>7r° — > 
K s K L n°. 

In our search for the process e + e~ — ► KsK^tt (|I5|) the decay channel K$ — > 2tt° was 

chosen. The following multi-photon reactions can be a source of background: e + e~ — > 
o 7r o 7r o 7 (g^ e+e - V1 _^ 37r o 7 ^ ? e+e - _^ ^3p ? e+e - _^ ^o^o _^ 3yr o 7 



W7T 

(0), e + e-^4 7 ,5 7 (0). 

The main selection cuts were E tot /2E > 0.38, iV 7 > 6, A^ cp = 0. The number of found 
7r°-s should be not less than 3 with two of them forming K$ meson, the recoil mass of 3 
7r°'s should be compatible with Kl mass. Fig.[59] presents the distribution over invariant 
mass for experiment and simulation. The contribution of QED background from processes 
(O) and from the process (Bl) was found to be negligible. To suppress possible background 



from the process flTB| ) with unknown cross section we used a set of harder cuts on ETON, 
PTRT, and 7r° triples with invariant mass close to that of rj meson. We estimated from 
VDM model, that under these cuts the expected contribution from this process should 
be small. The same hard cuts allow to suppress the process (|18f) |57| . Main background 
comes from the reaction e + e~ — > 07 with subsequent decay <fi — > KgK^. The number of 
such events was estimated [^3| to be (6 ± 3), which is close to observed 5 events. So, on 
the basis of the available data only upper limits of the cross sections in the energy range 
1300 1400 MeV can be placed: 

a(e + e~ -> 0vr°) < 0.2 nb (90% C.L.), 

a(e + e~ -> KK*) < 1.6 nb (90% C.L.), 

r(C -> e+e") • B(C -> 0tt o ) < 36 eV (90% C.L.), 



The obtained limits are lower than existing values [15 



7.5 Upper limit for electron width of /2(1270) meson 

C-even resonance production in e + e~ collisions is described in the lowest order in a by 
a Feinman diagram shown in (Fig. ^T]). The unitary limit of the electron width of the 
/2(1270) meson can be expressed through its two-photon width and a factor of ~ a 2 60fl , 
|5"I| ]. It can be estimated using the table value of r/ 2 ^ 77 



r u niUim.(/ 2 (1270) -> e + e-) ~ 3 ■ 10" 2 eV. (56) 

Actual cross section can be several times larger due to transition form factor. The only 
experimental result was obtained at VEPP-2M collider with ND detector ||. The e + e~ — > 
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/2(1270) — > 7r°7r° reaction ( |T9"D was studied in a 47 final state |)2], [15| and the following 
upper limit was obtained: r(/ 2 (1270) — > e + e _ ) < 1.7 eV at 90% confidence level. 

It this work, like in ||62|| , the reaction (|T9|) was studied again. The cross section was 



opproximated according to [63 



%{e+e- - / 2 (1270) - ttV) = 12.5 • ( • ^ " • *m 2 (20), (57) 

asZ \ m J (m — sj^ + rrrl z 

where s = 4.Eq, m, T, S ee B n o n o - are the /2(1270) meson mass, full width, branching 
ratios of the decays into e + e~ and 7r°7r° respectively. In the unitary limit the total cross 
section of the reaction (|i~9D (Fig. |62|) is about 1 pb at y/s = rrif 2 . This value is about 30 



% larger than that extracted from Breit-Wigner formula, used in [|62 |. 

Event selection For analysis of the reaction (|l~9|) , events with four photons in the 
final state where selected. The main background processes are the following: e + e~ — > 
47,57 (0) and e + e _ — > um° — > 7r°7r°7 (H). The cross section each of these processes is 
approximately three order of magnitude larger than that of the process (|T9"D. To choose 
selection criteria and estimate detection efficiency the Monte Carlo simulation was carried 
out for the processes (|T9"| , |T2"| , In order to suppress background the following selection 
criteria were applied: 

• Nj = 4; N cp = 0; 6 mm > 27°; E tot /2E > 0.85; P/E tot < 0.1; 

• C < 0; xlw < 15; 

• 85 MeV < 77114, m 23 < 185 MeV ; E 2 < 0.8E ; E 4 > 0.1E ; 

The shown above choise of certain combinations of photon pairs 777,14, ""7.23 is based on 
specific kinematics of the process ([19|) . The probability for other combinations to form 
7r°-s is about 20 %. The restrictions on the second and fourth photons energies - E 2 , E^ 
essentially reduced background from the process (|i~2"D. 

As a result of all described cuts, 3 experimental events left. The expected background is 
about 10 events mainly from the process @. To reduce background, the kinematic fit was 
applied for remaining events. The fit required the energy-momentum balance, presence of 
two 7T° meson in event. Events with {xio n o < 15) were selected. The experimental events 



are shown in Fig. |63|. The expected background from the processes(|T|, |J), estimated from 
simulation was found to be 1.0 ± 0.7 event. 

Results According to simulation the detection efficiency for the process ( |i~9"D is 21%. 
Based on the single found event the following upper limit could be placed: 

r(/ 2 (1270) - e+e-) < ■ h ■ T umtMm {f 2 {mO) - e+e~) = 0.4 eV, (58) 

* expected 

C.L. 90%, where k± = 3.89 is the Poisson coefficient for upper limit at 90% confidence level 
for one detected event. Obtained result in four times lower than previous one, obtained 
with ND [jTSQ, but still 15 times higher than unitary limit: r(/ 2 (1270) — > e + e~). 



8 Conclusion 

In Tables ^| and |7| the SND physical results, obtained up to now, are shown. Let us list 
the most important of them: 

1. The existence of meson radiative electric dipole decays — > / 7,a 7, announced in 
our previous work [0, is confirmed. The measured branching ratios and mass spectrum 
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Table 6: Measured branching ratios of particle decays in 0- 96 experiments. 



Decay mode 


Branching ratio 
(BR) (this work) 


Expected 
BR, ref. 


PDG,1996 or re- 
cent result 


Comments 


A\ v -77-0^77-0 r\i 

(JJ r /( /( V 


h 14 4- n in 4- 
0.12) • 10~ 4 


io- 4 ^io- 5 


< io- 3 


20| 


JTlloL UUocI Vd 

tion 


phi — > r]TT 0, y 


(8.3 ±2.3) • IO" 5 


lO^lO^ 4 




First observa- 
tion 


<p — > 7] 7 


(6.71^) • 10 5 ; 
< 1.1-10" 4 (in 7- 

y IIlULlt; j 


0.6 -10- 4 ^ 
10- 4 [18| 


(i-2±8:D • io- 4 
[0 


Second mea- 
surement 


cf) — > urn 


(5.7±fi) ■ 10- 5 


5- 

34 


IO" 5 J3^, 


< 5 ■ IO" 5 0, 


16| 


First observa- 
tion 


— > T]7T 7l 0/ J 


< 2 ■ 10~ 5 




L.5 
16 


-=-2-10~ 5 


< 2 ■ 10" 4 


First attempt 


(p — > ?77 


^l.Zl zt U.Uo It 

0.05)% 




(1.26 ± 0.06)% 
Q 


Most precise 
measurement 


— ► ?7e + e~ 


(1.42 ± 0.39 ± 
0.23) ■ 10" 4 


( 
i 

[ 


L.l 

fo 
or 
12 


• io- 4 

r unit 
m.f.) 


(1.3±85) 



(1.1±0. 
10~ 4 [14 


■ IO" 4 

5±0.2)- 
I 




— > e + e~7 


(6.8 ±1.1 ±0.7)- 
10~ 3 


6.3 • 10~ 3 


(4.9 ± 1.1) • 10~ 3 






6.9 • 10~ 5 


1.9 • 10~ 9 

m 


< 1.9 - IO" 5 




/a(1270) - 

e + e~ 


T < 0.4 eV 
(BR < 2.2-10- 9 ) 


r 

0.03 eV 
(Unit. limit) 


T < 1.7 eV 
|62 




C.L. 90%, the 
best upper 
limit 
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Table 7: The e+e annihilation processes under study in the energy range 
2E = 985 1400 MeV. 



The reaction 


The average 
cross-section 


Comments 


e+e~ — > c<J7r — > 7r + 7r~7r°7r° 


7.6 ±0.8 n& 


2E = 


e + e _ — > CJ7T — > 7T 7T 7 


0.65 ±0.04 ±0.04 nb 


2E = 


e+e - — > e + e~7 


23.6 ±0.3 nb 


2E = 985 4 1040 MeV; 9 min > 36°; > 5° 


e+e~ — > e + e _ 77 


0.43 ±0.09 nb 


2E Q = 985 4 1040 MeV; 9 min > 36°; Aip > 5° 




1.82 ±0.21 ±0.25 nb 


2£ = m ; min > 27°; Aip > 6° 


e+e~ — > 7r + 7r _ 7r + 7r _ 


3 + 24nb 


2£ = 1000 4 1380 MeV 


e+e - — > CJ7T — > 7r + 7r _ 7r°7r° 


3 4 15 nb 


2E = 985 4 1380 MeV 


e+e~ — > 7r + 7r~7r°7r° excl. un 


4- 15 nb 


2E = 985 4 1380 MeV 


e+e~ — > 7T + 7r _ 7r° 


2.4 4 4.4 n& 


2E = 1100 4 1380 MeV 


e+e" -> KsKl 


13 4 0.5 nb 


2E = 1040 4 1380 MeV 


e+e~ — > 07T° 


< 0.2 n& 


2E = 1040 4 1380 MeV; C.L. 90% 


e+e- -> ifX* 


< 1.6 nb 


2E = 1040 4 1380 MeV; C.L. 90% 



of 7T°7r and nn system support the 4-quark model structure of fo and ao mesons. 

2. The decay <fi — > 7/7 is observed, confirming the 1997 CMD-2 first observation of this 
decay, but SND branching ratio almost two times less, unless the errors is still high. 

3. The decay <fi — > un is observed for the first time. 

4. The measurements of n — > e+e~7 and — > ne + e~ decay branching ratios agree well 
with existing data. 

5. The cross-sections of e+e - — > 7r+7r~7r°, 7r+7r~7r°7r° n + n~ n + n~ reactions were measured 
in agreement with previous results. 

6. In analysis of e+e~ — > 7r+7r~7r° process the effect of p — ^-interference was observed in 
final state 7r+7r~ mass spectrum. 

7. We continued study of high order QED processes e+e - — > e+e~7, e+e-77, 777. 
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Figure 1: Week- by- week graph of integrated luminosity accumulated during experiments 
with SND. 




a - with wiggler, 1 985-86, ND 
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Figure 2: Averaged over live time luminosity of the VEPP-2M collider in ND experiments 
in 1983-^1987 and SND experiments in 1995-^1998. Solid line corresponds to the L ~ Eq 
law and normalized to L = 2 • 10 30 cm~ 2 s~ l at E = 500 MeV. 
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Figure 3: Integrated luminosity, accumulated with SND up to July 1998, as a function 
of energy. The histogram bin width is equal to 20 MeV. 




700 
m w , MeV 

Figure 4: Two-photon mass distribution 
in search for the process e + e" — * — > 
r]e + e~ , ?7 — > 77. Histogram - experimen- 
tal data; hatched histogram - simulation 
of the decay — > r]e + e~; dashed line - 
simulation of the QED process e + e~ — > 
e + e~77. 
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Figure 5: Optimal fit of the measured 
two-photon invariant mass distribution in 
search for the process e + e~ — > — > i]e + e~ , 

7/ 77. 
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Figure 6: Optimal fit of the energy de- 
pendence of e + e~ — > 4> — > ^e + e~, 77 — > 
77 visible cross section (for 500 MeV < 
m 77 < 600 MeV). 




Figure 8: The normalized energy spec- 
tra of the most energetic photon in experi- 
mental events of the reaction e + e~ — > <p — > 
777 —>■ 77 and simulated events of the re- 
action e + e~ — > (ft — > 77/7 — >■ 77. The cut at 
-Emaz = 0-Q5Eq is shown by arrow. 
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Figure 7: Optimal fit of energy depen- 
dence of <p — > 777 , 77 — » e + e~7 visible cross 
section. 
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T\%% mass, MeV 

Figure 9: The reconstructed mass spec- 
trum of 777r°7r° system for simulation of the 
process e + e~ — > <p —> ?77 - ► 37r°7 — > 77. 
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Figure 10: Two-dimensional plot of the 
measured mass of 77(550) meson versus its 
momentum in search for 777r°7r°7 events. 
Circles are simulated events of the reac- 
tion e + e~ — > — > 7/(958)7 — > 777r°7r°7, 
square represents one event from simula- 
tion of the process e + e~ — ► — > 777 — ► 
— > 77, triangles are experimental 
events. The vertical dashed line shows the 
momentum cut at i 3 ^ >50 MeV. 



Figure 11: Recoil mass distribution of 
the most energetic photon in search for the 
— > 7/7 decay. Circles with error bars - 
experimental data (11239 events); curve 
- optimal fit of the experimental data; 
histogram - estimated sum of contribu- 
tions from background processes e + e _ — > 

(f) -> TT + TT-TT (gl]), 
0^0 



7T 7T 7T 7T 



e e 



UJ7V 



5f); hatched histogram - con- 



tribution from the process e + e 
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Figure 12: Distributions of normalized 
energy deposition of neutral particles, in 
search for the decay — » rj'j. a) Cir- 
cles with error bars - experimental data; 
histograms - estimated contributions from 
background processes: hatched histogram 

- e + e~ — > (f> — ► 7T + 7r~7r° dotted line 

- e + e~ — > (f) — ► tun® — ► 7T + 7r-7r 7r° 
dashed line - e + e~ — > 177, ^ — ► 7r + 7r~7r° 
(|32j). b) Expected contribution from the 
process under study at Br((f> — > 7/7) = 
10" 4 . 



Figure 13: Distributions of normalized 
energy of the 2-nd photon in search for 
the decay — > 7/7. a) Circles with er- 
ror bars - experimental data; histograms - 
estimated contributions from background 
processes: hatched histogram - e + e~ — > 
(|3~ID; dotted line - e + e~ — > 



7T + 7T 7T° 



,0 



</> — > UJTX 

- e + e~ — > 777, 



7T + 7T 7T 7r° 



5j); dashed line 
tt+tt-tt (0). b) Ex- 
pected contribution from the process un- 
der study at Br(4> — > rfy) = 10~ 4 . 
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Figure 14: Distributions of experimental and simulated events over invariant masses 
of photon pairs in search for the — > 7/7 decay. Number of simulated events is not 
normalized. 
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Figure 15: Distributions of experimental (d) and simulated (a - c) events over recoil 
mass of the photon with minimal energy versus invariant mass of the pair of photons with 
maximal energy in search for the decay <p — > 7/7. Numbers of simulated events are not 
normalized. 
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Figure 16: Distribution of experimental 
and simulated events over invariant mass 
of two photons with maximal energy in 
search for the decay <p v'l- 
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Figure 18: Distribution of 7r°7 mass for 
m n7T < 700 MeV in search for — > 7r°7r°7 
decay. Points - data, histogram - sim- 
ulation, shaded histogram - sum of sim- 
ulated contributions from — > 777 and 
— > 7r°7r°7 decays, arrows - selection cuts 
for uo7T° class. 
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Figure 17: Distribution of experimen- 
tal and simulated events over recoil mass 
of the photon with minimal energy in 
search for the <fi — ► 7/7 decay. Hatched 
histogram and dashed line show the es- 
timated contribution of the background 
process e + e~ — ► 7/7, 77 — > 7r + 7r~7r (|32| ) 
and its linear approximation. Histogram 
and smooth curve show the distribution 
of experimental events and its optimal ap- 
proximation. 
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Figure 19: a, b - distributions of cosine of ^ - the angle between directions of tt° and 
recoil 7 in the rest frame of 7r°7r° system; c, d - distributions of 9 - angle of the recoil 7 
with respect to the beam. Points - data, histogram - simulation. 
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Figure 20: a - invariant mass distribution 
of 7r°7r° pairs for selected 7r°7r°7 events. 
Histogram - data, shaded histogram - 
estimated background contribution from 
e + e~ — > urn and <p — ► 777; b - detection 
efficiency for vr 7r°7 events. 
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Figure 21: Measured 7r°7r° invariant 
mass spectrum with efficiency corrections. 
Points - data, solid line - fit. 
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Figure 22: Energy dependence of the vis- 
ible e + e~ — > 7r°7r°7 cross section. Points 
- data, solid line - fit, dotted line - esti- 



mated background from e + e 
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Figure 23: Invariant masses of pairs 
of most energetic photons in search for 
cf) — > r/7r°7 decay: circles with error bars 
- experimental data, (a) - simulated sig- 
nal from <p — ► ?77r 7 decay corresponding 
to a branching ratio of 0.7- 10~ 4 , (b) - es- 
timated background from the e + e~ — > ujit 
and — > ^7, / 7 events, (c) - sum of (a) 
and (b). 
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Figure 24: Spectrum of the most ener- 
getic photon in an event in search for — > 
?77r 7 decay: a) photons quality £ < —4, 
b) < £ < 10. Circles with error bars de- 
pict experimental data, shaded histogram 
- simulation of the process (|5|) and clear 
histogram - simulated sum of (|5p and (0) 
with a Ei? = 10" 4 . 



Figure 25: Distribution of events over £ 



777r°7 decay: a) events 



in search for 
with E imax /E < 0.7, b) events with 0.7 < 
Erymax/Eo < 0.8. Circles with error bars 
- experimental data, shaded histogram - 
simulation of the process and clear 
histogram - simulated sum of (|5|) and (|2|) 
with a BR = 10" 4 . 
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Figure 26: Distribution of events over \ 2 a 
in search for — > r]7r°7 decay: a) events 
with £ < —4, b) events with < ( < 
10, circles with error bars - experimental 
data, shaded histogram depicts simulation 
of the process @ and clear histogram - 
simulated sum of (0) and (&) with a BR = 
10" 4 . 
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Figure 27: The recoil mass distribution 
of the most energetic photon in an event. 
Points - data, histogram - simulation, 
hatched histogram - simulation of the 
background decay 4> — ► K$Kl- 
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Figure 28: Distribution of 77 invariant mass closest to that of 7/-meson for 
photons not included into 7r° in search for the decay (f) — > r]7T j. a) his- 
togram - experimental data; hatched histogram - expected background; 
b) (f) — > ?77r 7 simulation. 
20 
Nl8 




17.5 20 

r\ny 

— > r]7i 0r y decay, with 

additional cut ( < —4. Points with error bars - experimental data; hatched histogram 
- expected background; histogram - sum of expected effect at £>(0 — > r/nj) = 10~ 4 and 
background contributions. 



Figure 29: Distribution of the XjL-o- parameter in search for 
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Figure 30: Distribution of the parameter 
( in search for the decay — > f]TT ^ with 
additional cut xL-o 7 < 7. Circles with er- 
ror bars - experimental data; hatched his- 
togram - expected background; histogram 
- sum of expected signal at B(<p — > 77717) = 
10~ 4 and background contributions. 
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Figure 32: Distribution of the ir/K sep- 
aration parameter in search for the pro- 
cess e + e~ — > urn — > 7r + 7r _ 7r°7r°. Arrow 
shows a cut suppressing background from 
cf) — ► K + K~ decay. 
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Figure 31: Distribution of the 7771"° invari- 
ant mass in a search for the — > 7/7r°7 de- 
cay. Solid line is a calculated distribution 
for the optimal value of fit parameters. 
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Figure 33: Distribution of the parameter 
X3tt for experimental and simulated events 
of the process e + e" — > ujtt — > 7r + 7r~7r°7r° 
(hatched histogram). A peak in the exper- 
imental data distribution is due to the the 
process e + e~ — > <p ~^ 7r + 7r _ 7r° ([31]) . Arrow 
shows the cut value for this parameter. 
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Figure 34: Distribution of the parame- 
ter m rec . for experimental data and sim- 
ulated process e + e~ — > ujti — > 7r + 7r~7r°7r° 
(hatched histogram). A peak at //-meson 
mass in experimental data is determined 
by contribution from the process e + e" — > 
— >■ 777, eta — > 7T + 7r _ 7r fl3"2|). Arrow 
shows the cut value for this parameter. 
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Figure 36: Distribution of the parameter 
Xatt for experimental (histogram) and sim- 
ulated (hatched histogram) events of the 
process e + e" — > uii c 
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Figure 35: Distribution of the parame- 
ter for experimental (histogram) and 
simulated (hatched histogram) events of 
the process e + e~ — ► 7r + 7r~7r 7r . Dots 
with error bars - simulation of the pro- 
cess e + e~ — > UJ7T° — > 7T + 7T~7T°7T° With 

7% contribution of the process e + e~ — > 
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Figure 37: Visible cross section of the pro 
cess e + e~ — > ujtt — > 7r + 7r~7r°7r° as a func- 
tion of energy and its optimal approxima 
tion. 
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Figure 38: Distribution of the parame- 
ter X7r°7r°7 f° r experimental (histogram) 
and simulated events of the processes 
e + e~ — > un — > 7r°7r°7 (hatched his- 
togram), e + e~ — > — > 777 — > 37r°7 
(dashed-line histogram). 
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Figure 40: Distribution on the param- 
eter x-ymax for experimental (histogram) 
and simulated events of the processes 
e + e~ — > cj7r° — > 7r°7r°7 (hatched his- 
togram), e + e~ — >■ — > ?77 
(dashed-line histogram). 
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Figure 39: Distribution of the param- 
eter for experimental (histogram) 
and simulated events of the processes 
e + e~ — > urn — > 7r 7r°7 (hatched his- 



togram), e + e — > - 
(dashed-line histogram). 
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Figure 41: Visible cross section of the pro- 
cess e + e~ — > CJ7T — ► 7r°7r°7. The optimal 
fits of visible cross section and resonance 
background are shown. 
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Figure 42: Comparison of experimental data and simulation for the process e + e~ — > 
e + e~7: a) normalized energy of the charged particle; b) normalized photon energy; c) 
angle between charged particles; d) angle between photon and nearest charged particle. 
Triangles - experimental data, histogram - simulation. 
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Figure 43: Comparison of experimental data and simulation of the process e + e~ — > 
e + e~77: a) normalized charged particle energy; b) normalized photon energy; c) angle 
between charged particles; d) angle between photons; e) minimal angle between photon 
and charged particle; f) invariant mass of photon pair. Circles - experimental data, 
histogram - simulation. 
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Figure 44: Visible cross section of the pro- 
cess e + e~ — > 777. Upper curve - optimal 
fit. Lower curve - resonance background. 
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Figure 46: Normalized spectrum of invari- 
ant mass of photon pairs in the process 
e + e~ — > 777, three entries per event. Cir- 
cles with error bars - experimental data; 
histogram - simulation. 
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Figure 45: Distribution of normalized 
minimal photon energy in the process 
e + e~ — > 777. Circles with error bars - ex- 
perimental data; histogram - simulation. 
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Figure 47: Distribution of azimuth 
acollinearity angle between 1-st and 2-nd 
photons in the process e + e~ — > 777. Cir- 
cles with error bars - experimental data; 
histogram - simulation. 
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Figure 48: Experimental cross section of the process e + e — > it + it it + it . Only statistical 
errors are shown. 
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Figure 49: Distribution of invariant mass of 7r + 7r~ system in the process e + e~ - 
7r + 7r~7r + 7r~. Circles with error bars - experimental data; histogram - LIPS simulation. 
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Figure 50: Distribution of 7T°-meson recoil mass in the process e + e~ — > 7r + 7r~7r°7r°. Circles 
with error bars - experimental data; hatched histogram - LIPS simulation; histogram - 
simulation of the process e + e~ — > u;7r — > 7r + 7r~7r°7r°. 
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Figure 51: Measured cross section of the process e + e~ — > 7r + 7r~7r°7r° . 
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Figure 52: Measured cross section of the process e + e~ — > 7r + 7r~7r°7r°. 
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Figure 53: Distribution of invariant mass 
of 7T° pairs in study of the process e + e _ — > 
K$Kl. Points with error bars - experi- 
mental data; histogram - simulation. 
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Figure 54: Cross section of the process 
e + e~ — > KsKl: • - SND experimental 
data, o - OLYA experimental data, trian- 
gles — calculation. 
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Figure 55: Energy dependence of the to- 
tal cross section of the process e + e _ — > 
tv + tv~tv q ( |5DD measured with ND detector. 
Smooth curve - vector dominance model. 
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Figure 57: Comparison of experimental 
results obtained with the ND (clear tri- 
angles) and SND (black triangles) detec- 
tors for the cross section of the process 
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Figure 56: Distribution of for se- 
lected experimental events of the process 
e + e~ — > tv + tv~tv°. Arrow shows the cut-off 
for this parameter. 
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Figure 58: Observation of the p-oj in- 
terference in invariant mass spectrum of 
charged pion pairs in the reaction e + e~ — > 
tv + tv~tv°. Histogram - experimental dis- 
tribution of tv + tv~ invariant mass; hatched 
histogram - experimental distribution of 
7r°7r ± invariant mass; black histogram - 
simulated tv + tv~ invariant mass distribu- 
tion for the background process e + e~ — > 
tv + tv~tv°tv . 
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Figure 59: Distribtion of selected events over recoil mass of the 7r°7r 7r° system versus 
invariant mass of 7r 7r pair in the search for e + e~ — > K s K l t: decay. 
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Figure 60: Distribution of experimental 
events over E to t/2E versus P/E tot , in the 
search for Kg — > 3tt° decay. 
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Figure 62: Calculated total cross section 
of the process e + e~ — ► / 2 (1270) — ► 7r°7r°. 




Figure 61: Feinman graph for direct pro- 
duction of C-even resonance in e + e~ colli- 



sion. 
> 



180 
160 
140 
120 
100 















• 




* 










* 
* 














" I 


, , , 


, , , 


, , , 


, , , 



100 120 140 160 180 



m 1 4 , MeV 



Figure 63: Scatter plot of invariant masses 
of photon pairs in search for the process 
e+e" -> / 2 (1270) -> Small dots - 



simulation of the process (PUf); triangles - 
simulation of the background process dl2|) ; 
stars - simulation of the background pro- 
cess (^|); circle - the only experimental 
event, which passed all the selection cuts. 
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